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I , INTRODUCTION 


1,1 Oblectives and General Relevance of the Study 

Three dimensional turbulent boundary layers are experienced in a 
wide variety of flows over both stationary and rotating bodies, such as 
the blades of fluid machinery. The crossflows, even though in some cases 
small, significantly affect the boundary layer growth and its structure. 
In cases where the flow undergoes large turning and has crosswise 
pressure gradients (e.g. on the leading or the trailing edge of a stator 
blade) or experiences in addition rotation effects (e.g. rotor blades), 
the cross flows are large and any two-dimensional simplification of the 
flow will be grossly in error. 

In addition to producing cross flows, the rotation and curvature of 
a body alter the basic turbulence structure of the flow. For example on 
the convex surface of a bend, the flow tends to stabilize and even 
re-laminarize for large values of curvature. On the other hand, the 
flow is destabilized on a concave surface and may even separate for 
large values of curvature. Similarly rotation also stabilizes or 
destabilizes the flow, depending upon its direction with respect to the 
flow. 

Lack of knowledge on these three-dimensional turbulent boundary 
layers, especially those with rotation and curvature, has hampered the 
technological progress in areas such as the fluid machinery (e.g. 
compressors, pumps, propellers, turbines, etc.). The rotation induced 
Coriolis and centrifugal forces generate spanwise flows inside the 
blade boundary layer, the magnitude of which depends on the angular 
velocity of the blade, the flow coefficient, the rotor blade geometry 




and Its cascade parameters such as space/chord ratio, aspect ratio, 
stagger angle, etc. This spanwise migration of the flow results in a 
decrease of the boundary layer growth near the hub and an increase near 
the tip, an increase in skin friction and changes in the turbulent 
characteristics of the flow. 

The flow losses due to blade boundary layers form a major portion 
of the total aerodynamic losses in the turboraachinery stage. It varies 
directly with free stream velocity, boundary layer thickness, skin 
friction coefficient, magnitude of cross flow inside the boundary layer 
and location of flow separation point (if any) . All these variables in 
turn are strongly dependent on the three-dimensionality of the flow. 
Hence omission of the three-dimensional effects due to rotation and 
curvature in the prediction of aerodynamic efficiency on the basis of 
stationary two-dimensional boundary layer theories and cascade data is 
entirely inadequate. 

In addition, near the tip, complex flow interactions take place 
between the blade boundary layer and the annulus wall. The outward 
spanwise flow in the blade boundary layer is deflected inward by the 
annulus wall to the center of the passage. This results in an inward 
spanwise flow with a marked reduction of mainstream velocity. The shape 
of velocity profiles in this "interference layer" (so named, because it 
is formed by the interference between the annulus wall boundary layer 
and the blade boundary layer flow near the tip) is unconventional and 
resembles a wake. The radial inward flow component Increases and the 
mainstream flow component decreases with Increase in distance from the 
blade. The direction of turbulent shear stresses is likewise found to 
reverse in accordance with the mean velocity gradients, and the turbulent 
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intensity level also increases with distance from the blade. This 
results in an increase in the mixing losses and the turbulence generated 
broad band noise in rotating fluid machinery. It also explains some of 
the discrepancy in the end wall loss correlations between the measure- 
ments in a rotor and correlations based on stationary cascade data. 

It is thus obvious that flow analysis and measurements in an actual 
rotating fluid machinery are essential, both in studying the effects of 
rotation and curvature on three-dimensional turbulent boundary layers 
with arbitrary pressure gradients, and in incorporating their effects 
in the design of turbomachinery blading. An extensive and systematic 
investigation was therefore conducted to study these important and 
basic differences between stationary and rotating turbulent boundary 
layers on the blades of an axial flow fluid machinery rotors, and in 
particular, on the blades of an axial flow rocket pump inducer. An 
inducer has long and narrow blade passages and therefore represents the 
worst case of three-dimensional effects. The previous study was 
performed on a single helical blade of large chord length, rotating in a 
fluid, otherwise at rest, and enclosed in an annulus. This represents 
a study of three-dimensional boundary layers with zero pressure gradients 
and includes the effects of rotation and curvature only. The interference 
due to adjacent blade boundary layer and the annulus wall in the ducted 
rotating passages formed by the rotor blades in actual fluid machinery 
could not be simulated. The present study was conducted inside the 
straight helical channel formed by the blades of a four bladed flat 
plate axial flow inducer. Even though the direct effects of cambered 
blades of actual turbomachinery blading could not be investigated 
experimentally, the theoretical analysis presented here is valid for 
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axial flow fluid machinery rotors with arbitrary blade geometry 
and blade loading. 

The following objectives were kept in mind in conducting the study: 

1. To study the effects of rotation and curvature on three- 
dimensional turbulent boundary layer characteristics such as mean 
velocity profiles, skin friction, turbulence structure, boundary layer 
growth, etc. in a rotating channel with arbitrary pressure gradients. 

2. To study the Interference effects due to the presence of an 
annulus wall or adjacent blade on rotating blade boundary layers in an 
actual fluid machinery. 

3. To develop calculation methods (e.g. integral method) to 

predict the gross three-dimensional boundary layer parameters 6, , , £ 

11 w 

and H. These are essential for calculating the boundary layer flow 
losses and in predicting the overall aerodynamic performance of fluid 
machinery. 

4. To develop probes, probe traversing device, signal processing 
system and Instrumentation for accurate measurement of boundary layer 
flow characteristics in a rotor frame of reference. 

1.2 Statement of the Problem 

This is an experimental and analytical study of the effects of 
rotation and curvature on the three-dimensional turbulent boundary layer 
flow characteristics in a rotating channel with arbitrary pressure 
gradients. The following characteristics of the boundary layers are 
studied throughout the rotating channel except where otherwise stated: 

1. Mean velocity profiles 

2. Turbulent intensity and shear stress components 
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3. Spectral distribution of all six turbulent correlations of 
Reynolds stress tensor (at mid-chord and near the blade tip only) 

4. Growth of boundary layer momentum thickness 

5. Limiting streamline angle or skew angle 

6. Shape parameter 

7. Skin friction and its correlation with other boundary layer 

parameters (0- - » e and H) . 

11 w 

In addition, the interference effects on the boundary layers due 
to the presence of end walls and boundary layers on the adjacent blades 
are also studied. 

1. 3 Previous Related Work 

Exact solutions of Navier-Stokes equations for rotating bodies have 

been obtained only for the simplest cases of three-dimensional rotating 

laminar boundary layers, Karman (1) obtained the complete solution for 

laminar flow over a rotating infinite disk immersed in an infinite 

stationary fluid using similarity analysis for the velocity profiles. 

Bama (2) and Benton (3) have obtained the laminar solution in rotating 

pipes by assuming the flow to be a small perturbation of Hagen- 

Poiseuille flow, and found that the friction factor "f" is increased by 

rotation given by f/f^_Q= f (Re^ ,/Ro"^) , where Ro^ is a rotation parameter 

=* i^d/U^ and Re^^ is Reynolds number based on hydraulic diameter and 

centerline velocity U^, d is the pipe diameter and ^2 is the angular 

velocity of rotation. Experimental results of Treffethen (4) in the 
2 3 

range of 2 x 10 < Re, , Ro, < 3 x 10 confirmed the above conclusion. 

h d 

The stability analysis of Chawla (5) for laminar rotating boundary 
layers has shown that the Tollmeln-Schlichting waves are not affected by 
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rotation but Taylor-Gortler vortex cells are damped and the transition 
to turbulence is delayed. 

Horlock and Wordsworth (6) and Banks and Gadd (7) Investigated the 
three-dimensional laminar boundary layer characteristics on a turbo- 
machinery blade. Similar Investigations on a helicopter blade have been 
reported by Fogarty (8) and McCroskey and Yaggy (9). Effects of 
curvature and pressure gradients are neglected and small cross flow 
assumptions are made. They conclude that the rotation decreases the 
boimdary layer growth compared to the stationary case, primarily due to 
spanwise migration of the flow. 

For turbulent boundary layers on rotating bodies, lack of information 
on the exact nature of turbulent shear stress terms in the equations of 
mean motion (known as the closure problem in turbulent flows) has eluded 
an exact solution even in the simplest case of incompressible flow over a 
rotating disk. Approximate solutions based on momentum Integral techni- 
que have been obtained for a few cases but the empiricism involved is 
based on scarce experimental data in rotating systems. The procedure 
involves the integration of the differential equations of motion across 
the boundary layer and thus eliminates the need for the nature of 
turbulent stresses which vanish both at the wall and outside the 
boundary layer. The final equations contain momentum thickness 
parameters, which are related to each other and solved by assuming 
similarity in velocity profiles and an empirical relation for the skin 
friction. The following studies of this nature have been made: 

Turbulent Boundary Layers on Axisymmetrlc Bodies . Karman (1) was 

first to obtain a solution for turbulent boundary layer growth on a rota- 
ting disk. He used l/7th profile for the mainstream velocity, skin friction 
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coefficient based on stationary flat plate data and a constant value 
of 0^, the limiting streamline parameter. The parameter Is defined 
as the tangent of the angle which the local streamline makes with the 
Invlscld free streamline, as the distance of the former from the blade 
surface approaches zero. Cham and Head (10) obtained a numerical 
solution of the momentum Integral equations for turbulent boundary layer 
growth on an axlsymmetrlc disk and concluded that Is nearly Invariant 
with radius. The atmospheric turbulent boundary layer (11) also falls 
In this category. Solutions have been obtained (e.g. Blackadar and 
Tennekes (12)) by assuming that the Rossby Number Is very large. The 
Rossby Number, Ros, Is the Inverse of rotation number Ro and for 

atmospheric boxmdary layers. Its value Is Indeed very large compared to 

1 

one. However, In most of the flows In laboratory and Industrial 
rotating equipment, the Rossby number Is small and the large Rossby 
number assumption presents a serious limitation on the practical 
application of this theory. 

Turbulent Boundary Layers in Centrifugal Impellers . The axis of 

rotation in centrifugal impellers is perpendicular to the mainstream 

flow direction and hence large Coriolis forces are produced normal to 

the plane containing the axis of rotation and the flow. Hill and Moon 

(13), Moon (14), Moore (15), Haleen and Johnston (16) have investigated 

such flows. Moon's (14) result for rectangular channel rotating about 

an axis perpendicular to the flow showed that the boundary layer growth 

1 

on the low pressure side is increased substantially and the velocity 
profiles are of power law type, while on the high pressure side the 
boundary layer is much thinner and the velocity profiles fit the 
universal law of the wall. Skin friction measurements showed a decrease 
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on the suction side and substantial increase on the pressure side. 

Moore (15) studied the flow on the center line of the rotating ducts of 

various aspect ratios (ratio of the two sides of the rectangular section) 

ranging from 0,5 to 7.33. His results show an increase in friction 

factor of *^207o on the leading side and a similar decrease on the trailing 

side. However uncertainty in data is too high to allow generalization to 

predict the effects of rotation in such flows. Haleen and Johnston (16) 

did a detailed experimental study on the structure of turbulent boundary 

layers in a rotating rectangular duct. Using hydrogen bubble technique, 

they showed that on the suction side, a tendency toward supression of 

turbulence exists and the tendency increases with rotation number and 

for very large Reynolds number, the turbulent flow may even re -laminar ize. 

The rotation number, Roo, is defined here as ROp = 2f26/U_ where Q is the 

0 6 =0 

magnitude of angular velocity of rotation, 6 is the boundary thickness 
and is the free stream velocity. The rotation number is the inverse 
of Rossby number (Ros = U/2f^L). The former is more convenient to use in 
turbomachinery applications, while the latter is well known for its 
application in the geostrophic flows. 

Turbulent Boundary Layers in Axial Flow Machines . To date, very 
little work has been done to investigate the characteristics of rotating 
turbulent boundary layers on the blades of an axial flow fluid machinery. 
In many cases where the blade chord and the stagger angle are small, 
the main stream flow direction is almost parallel to the direction of 
rotation. The Coriolis force is small and very little secondary flows 
are developed due to rotation. However if the blade chord and the 
stagger angle are large, the boundary layers are thick. The secondary 
flows are large and the effect of Coriolis force cannot be neglected. 
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Banks and Gadd (7) investigated the turbulent boundary layer growth 
near the leading edge of a single helical blade of large stagger 
rotating in a stationary fluid. They approximate the blade by a flat 
sector of a disk with leading and trailing edges. The two momentum 
integral equations were solved by assuming a series solution, small 
cross flow and l/7th profile for the mainstream velocity. Later, 
Lakshminarayana et al. (17) conducted an extensive experimental and 
analytical investigation of the three-dimensional turbulent boundary 
layer growth on a single rotating helical blade of a rocket pump Inducer. 
Using the momentum Integral approach, they predicted the boundary layer 
growth and the limiting stream line angle of the flow over the entire 
blade. They also obtained the nature of mean velocity profiles in the 
outer part of the turbulent boundary layer and a skin friction 
correlation to show the direct effects of rotation. The skin friction 
stress was found to Increase with rotation according to the formula 

c /c = 1 + c (e X /r) 

where is the limiting streamline angle and x is the distance from 
the leading edge. 

Effect of Curvature in Boundary Layers . The effect of curvature on 
boundary layers has been discussed by Bradshaw (18), Stratford et al. 
(19) . Stratford et al. proposed an expression for the change in mixing 
length, based on classical mixing length arguments, whereas Bradshaw 
drew an analogy between streamline curvature and buoyancy. A first order 
effect of curvature, suggested by Bradshaw on the basis of Monin -Oboukhov 
meteorological formula is 


1 /2 

where il = /(S'i/Sy)- Ri is Richardson Number = 2s (1 + s) 

where s = (u/R/3u/8y) and 3 i« a constant. Bradshaw suggests a value of 
3 =» 7 for the convex surfaces and ^ ^ for the concave surfaces . 

Bradshaw applied the same analogy to turbulent flows with rotation and 
concluded that the equivalent Richardson number for the case of rotation 
can be written by replacing (u/R) with 2f^y, i.e. 

Ri = “■ (2Qy) /(9u/3y) 

where f^y is the component of rotation vector ^ normal to the plane of 
mean rate of strain. Experimental data of Haleen and Johnston (16) for 
the rotating case agree fairly well with Bradshaw’s predictions. The 
combined effects of rotation and curvature on turbulent quantities can 
be predicted by using a combination Richardson number. 

The experimental data of various authors (14, 15, 16, 17, 20) in 
the rotor blade passages of centrifugal and axial-flow turbomachinery 
has clearly shown that rotation and curvature greatly modify the boundary 
layer growth, turbulence and skin friction characteristics. But a 
generalized analytical and experimental study to explain these changes 
in flow behaviour, especially in mean velocity profiles and turbulence 
quantities which are so important in predicting the aerodynamic and 
acoustical performance of a turbomachinery rotor, was lacking and 
prompted this investigation. 

1.4 Method and Means of Investigation 

As discussed in previous sections, three-dimensional flow analysis 
and measurements in an actual fluid machinery are essential to study 
the combined effects of rotation and curvature on the characteristics of 
three-dimensional turbulent blade boundary layers with arbitrary pressure 
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gradients. Analytical expressions for mean velocity profiles are 
derived in Chapter II for the inner, outer and interference regions of 
the boundary layer in generalized orthogonal curvilinear rotating 
coordinate system suitable for an axial flow turbomachinery rotor with 
arbitrary blade geometry, cascade parameters and blade loading. Boundary 
layer growth, limiting streamline angle and shape factor are predicted 
throughout the rotor channel using the momentum integral technique and 
by solving the resulting three first order partial differential equations 
simulatneously (two momentum integral equations and an auxiliary 
entrainment equation ) by a finite difference scheme. 

The experimental work is carried out in a four b laded flat plate 

5 

inducer at a Reynolds number of 6.67 x 10 (based on peripheral velocity 
and radius at the tip) , which assures that the flow is turbulent over 
most portions Of the blade channel. Experimental test facility, data 
transmission and instrumentation are described in Chapter III. Extensive 
three-dimensional flow measurements of the blade static pressure, 
limiting streamline angle, mean velocity profiles, turbulence intensities, 
turbulent shear stresses, skin friction have been made throughout the 
rotor channel. Data reduction procedures and methods of measurements are 
reported in Chapter IV, 

Measurements were performed both at open throttle and partially 
throttled conditions. The former gave mild chordwise pressure gradients 
so that the major effects of rotation and curvature could be discerned. 
With throttling, large pressure gradients are produced with flow condi- 
tions comparable to an actual fully loaded axial flow rocket pump inducer 
studied by Lakshminarayana (21) . Laminar to turbulent transition and 
spectral measurements of all the turbulent quantities at mid chord and 
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near the tip are reported. Performance characteristics and exit flow 
measurements have also been conducted and the experimental results of 
all measurements are reported In Chapter V. They are compared with 
theoretical predictions of Chapter XI. Subsequent chapters Include 
discussion, conclusions and recommendations for future work and methods 
of Incorporating these effects in design of turbomachine blading. 
Direction cosines and dimensions of sensors, supports, etc. of the 
special boundary layer triple sensor hotwire probe, response 
characteristics of two dual channel hotwire anemometers and turbulent 
scales of the flow ect. are discussed in Appendix A. A systematic 
analysis of errors in measurements, their elimination and the accuracy 
of data Is given in Appendix B. The equations of Reynolds stresses In 
reference rotating coordinate system (x, y, r) and the qualitative 
effects of rotation and curvature are reported in Appendix C. 
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II. THEORETICAL ANALYSIS 


2,1 Qualitative Nature of the Flow 

Extensive flow visualization Inside the blade passages of an axial 
flow Inducer, reported by Lakshmlnarayana (22), show that the qualitative 
nature of the flow in rotating passages is different from their 
statlotiary counterpart. Inside the boundary layer of a rotating blade, 
there Is an Imbalance between the radial pressure gradient and the 
centrifugal force. This produces radial outward or inward flow 
depending on whether the latter force in a fluid machine is larger than 
the former or vice versa. 

In pumps, compressors, etc., (Figure 1) the absolute and relative 
tangential velocities are in opposite direction. The centrifugal force 
is therefore larger than the radial pressure gradient inside the blade 
boundary layer and hence outward radial flows are produced. On some 
locations of a turbine blade or stator blade, the centrifugal force is 
smaller than the absolute tangential velocity component and hence 
radial inward flows are produced. The inducer, which has the absolute 
and relative tangential velocities in opposite direction represents the 
first case, consequently, the radial flow in its blade boundary layers 
is outward. This spanwise migration of the flow, produces thinner 
boundary layers near the hub and thicker boundary layers near the tip. 
Since the blade passages of an inducer are long and narrow (chord/ 
spacing - 3.6), the inducer represents the worst case of the three- 
dimensional boundary layer. The cross flows are large and limiting 
streamline layers are thick and merge with each other after nddchord 
at most of the radial locations. Near the tip, due to the interference 
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YTjjt Radial Inward 
Flow Region 



Prol 



1. The Inner Layer *- y < 10, where the viscous forces dominate 
and the inertia forces can be neglected. 

2. The Buffer Layer/Inertial Sublayer - 10 < y <10 , where all 
the forces including Inertia, viscous and turbulent forces are important. 
A knowledge of the nature of turbulent quantities, especially the 
turbulent stresses is required to predict even the mean velocity profiles 
of the flow field. 

3 

3. Outer Layer - y > 10 , where viscous and turbulent stress 
gradients may be neglected. 

4. Interference Layer - It is formed near the tip in the central 
region of the channel section due to the interaction of blade boundary 
layers with the annulus wall. The outward radial flow near the blade 
surfaces is deflected inward by the annulus wall, to the center of the 
channel section. The resulting inward radial flow interacts with the 
malnflow, causing a reduction in streamwise velocity. Thus the shape of 
the mean velocity profiles is reversed in this region compared to those 
near the blade surface. 

2. 2 Flow Model and Governing Equations 

The flow of a Newtonian fluid is governed by the Navier-Stokes 
equations. In the boundairy layer on a rotating blade of a turbomachinery, 
the flow is inherently unsteady when viewed from a stationary coordinate 
system. It is steady, in the absence of inlet distortions, when viewed 
from a non-inertial coordinate system rotating at the same angular 

velocity as the rotor. The rotation introduces two additional 
acceleration terms viz, Coriolis acceleration (2SxW) normal to both the 
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relative flow velocity vector (W) and the rotation vector («) , and 
centrifugal acceleration (V(-“ (j?^xr) • (Sxr) . 

The Navier-Stokes equation in rotating frame of reference is, 

2fixW = - Vp*/p-vVx(VxW) (1) 

where p* = p + pgh-(p/2) (^xr) ,(S^r) , is known as the reduced pressure and 
W is total relative flow velocity vector, g is the acceleration due to 
gravity, and h is the height of the fluid element from a reference axis. 

The Coriolis and centrifugal forces modify the local force balance 
between inertia, pressure and viscous forces on a fluid element and are 
responsible for the observed changes in the mean velocity profiles, 
boundary layer growth, turbulent structure, skin friction stresses, etc., 
compared to an equivalent stationary flow system. 

The generalized equations given below are based on the following 
assumptions. These equations are not necessarily restricted to inducer 
blades and are valid for any axial-flow turbomachine with arbitrary 
pressure gradients . 

(a) The flow is incompressible, homogeneous, single phase and 
steady in a relative frame. 

(b) The blade elements are nearly radial; leading to the choice 

of curvilinear orthogonal coordinate system (x, y, r) shown in Figure 3. 

X is along the blade in the direction of the Invlscid streamline and y 
is normal to the blade. Both x and y are curved and orthogonal to each 
other as well as the radial coordinate r. (i.e. x and y lie in the 
co-axial cylindrical plane of the machinery). This system is related to 
cylindrical coordinate system through the equations. 
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dx = rd6 sin3 + dz cos3 


( 2 ) 


dy * -rd0 cos 3 + dz siti3 

(c) Only the dominant shear stresses are retained and the usual 
boundary layer approximations are made. 

Hence the principal curvatures encountered in this flow system are 


K 

xr 


sln^3 

r 


K 

yr 


K 

xy 


3x 




yx 3y 


(3) 


where, for example, is the principal curvature on the surface r =» 
constant in the direction of the parameter x. 3 is the streamline angle 
(x direction) measured from the axial direction (z). With the usual 
boundary layer approximations, the time averaged equation of mean motion 
in a rotating (x, y, r) system can be written as (Figure 3) , 


X momentum: 




+ w 


3u 

3r 


+ sin^3 

r 


w sin3 


+ 



1 ^ ^ I 
p 9x P 9y 


r momentum: 


9w , 9w . 9w . 2 q 

u-;^ + v-r-+w-r sin 3 

9x dy dr 
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~ + 2^2 u sin3 
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- 2f2 V cos3 = - 


1 9p* 
p 9r 
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(5) 


continuity: 
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where p*, the reduced pressure, neglecting the term due to gravity 
(pgh) , is 

P* = p - 'I' P (f2xr) • (S^xr) 

u, V, w are mean velocities in x, y, and r directions, ^ is the angular 

velocity, and are total viscous and turbulent shear stresses in 

X and r directions respectively. For small cambered blades or the flat 

plate inducer used in the experimental program, retaining the principal 

2 

curvature term = sin 3/r equations (4) , (5) , and (6) are simplified 


to: 


... . 3U - 9t 

, du , ow • o • 2. uw e 1 X 

u h V 1 w 7 ^ 2s2w sin p + sin p U — 

9x 9y 9r r e dx p 9y 


(7) 


2 

3w . 9w . 9w . 2^ sin B s ,0 

u -r H V 7T“ + w 7r~ + (U - u ) — 2J2(U - u) sinB = 

dx dy dr e r e 


. 9x 

^ 

p 3y 


( 8 ) 


9u , 9u , 9w , w 
F ^ r 


0 


(9) 


The pressure gradient term (9p /p9x) is replaced by its value in 

the external stream (U^ 9U^/9x) . Similarly, the term 9p /p9r is 
2 2 

replaced by (U^ sin ^/v - sinB) • The velocity and its gra- 

dient has to be derived from a three dimensional inviscid program 
(e.g. Cooper and Bosch (24)). 

A complete solution of these highly nonlinear equations has not been 
tried due to lack of knowledge of the nature of turbulent stresses in 
the rotating three-dimensional turbulent boundary layers. However, a 
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possible approach using Bradshaw's three-dimensional boundary layer 
program (25) Is outlined In Chapter VI and explained by Anand et al. 
(20) . It Is hoped that the present experimental data and that of 

[r- 

Gorton (23) on turbulent stresses would provide a basis to check this 
approach and supply the required Information for the complete solution 
of the equations of motion. 

In the following section, a solution for the nature of mean 
velocity profiles In the Inner, outer, and Interference layers of a 
rotating turbulent boundary layer Is derived from the equations (7) to 
(9) using the previous qualitative and quantitative Information on such 
boundary layers. The buffer layer, which requires a complete solution 
of equations of motion Is not tried due to lack of knowledge of the 
nature of turbulent stresses In this region. 

2.3 Velocity Profile Relationship 

Since the cross flow velocity component (w) Is coupled to the main 

stream velocity component (u) through the equations of motion. It Is 

anticipated that similarity in velocity profiles should exist if the 

velocities are normalized by the local velocity scale, and distances by 

the proper length scale. The two convenient local parameters of the 

flow are T (the shear stress at the wall) and U (the local mainstream 
0 e 

velocity at the edge of the boundary layer) . Using the properties of 

the fluid p and y, a length scale (v/u^) is used, where v is the 

kinematic viscosity and u. is the shear velocity defined as u. = A /p . 

” ”0 

For the case of Inner layer (being closest to the surface) , the proper 
normalizing velocity scale is the shear velocity u^ and for the outer 
layer (nearest to the inviscid free stream region, Figure 2), it is the 
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I 


local free stream velocity U^. If similarity exists, the normalized 
velocity components would be unique fimctions of normalized distance 
y"** = u^y/v from the surface. Analytical relationship between the two 
velocity components in these regions follows: 


2-3.1 Inner Layer (Region I - Figure 2) 

+ 

In the inner layer (I) , close to the wall (y < 10) viscoiis forces 
dominate and due to the no slip condition at the wall, the Inertia terms 
in the momentum equations (7) and (8) become very small and are therefore 
neglected. Therefore in this region, the Coriolis force and the 
pressure gradient balance the shear stress gradient. Equations (7) and 
(8) are simplified to 

-2nwslnS-U^^=i^ (10) 

20 u sine - ^ (^2^ - sln^e) = ^ (U) 


where k * U^/fir. If we assume that the stresses are completely l ami nar, 
the stresses can be replaced by mean velocity gradients. This is true 
in the case of a laminar boundary layer and the laminar sublayer (y"*” ^ 5) 
of a turbulent boundary layer. However, contribution due to gradient of 
turbulent shear stresses is negligible in this region. The logarithmic 
law of the wall which is experimentally shown to be valid for y"*" < 200, 
gives a zero contribution to the shear stress gradients in this region, 
if the estimate is based on the mixing length hypothesis, given by 


with 


■ ky and u/u^ 


1/K^ln y"*" + C 


(13) 
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3t /3y = 0. Here the subscript t denotes turbulent. Similarly 
*t 

3t /3y = 0. The proper velocity and length scales in this region are 
u^ and v/u^. Assuming similarity in velocity profiles, w and u may be 
written as: 

w/u^ = f(y^)* (14) 

= gCy'*') 

where y = Reynolds number based on shear velocity. 

Replacing the stresses by corresponding velocity gradients and 
using Equation (14), Equations (10) and (11) become, 

- 2 f - TG/r’ = g" (15) 

2 0.^ g - (2sin 3/k - sin^$)/r* - f" (16) 

where = f2vsin6/u^^, r’ = r^ “ (r/U^) 3U^/3x. At the 

wall, shear stress gradient must be completely balanced by the pressure 
gradient. This and the no slip condition results in the following 
boundary conditions on the blade surface. 

1. f(o) - g(o) = 0 

2. f’^(o) + g»^(o) - 1 

3. f"(o) = - (2sin3/k - sin^3)/i^* 

4. g"(o) = - TG/r’ 

The solution of f and g is obtained in a series form by writing 

2 3 2 

f = + a^ y + a^ y + a^ y +....+ b^/y + b 2 /y + . . . 

(17) 

2 3 2 

g = y"^ + 5^2 + 5-3 y^ +«...+ m^/y"^ + m2/y^ + . . . 

(18) 
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The various constants In Equations (17) and (18) are evaluated from the 

boundary conditions by substituting for f and g. This gives, 

a=b-=b_»m,=m_*»£.=0 
o 1 2 1 2 o 

f • (o) = - e„/ / 1+e^^ 

f"(o) = 2a^ = - (2sln /k 

g”(o) = 2^2 = -TG/r' 

where e = tan a , a » limiting streamline angle as defined In Figure 3 
w w w 

By successive differentiation of Equations (15) and (16) , it is 
found that the nth derivative Is related to lower derivatives as follows 

g = (-2J2 )f = -(2^2 ) g 


(19) 


-sin 0)/r* 


(n >2) (n > 4) 


( 20 ) 


and 

^n ^"2 ,*^o.2-n-4 

f = (+2f^ )g = -(2S2 ) f 


(n >2) (n > 4) 

All the coefficients of f and g can be determined by taking the nth 
derivative at the wall and relating it to known coefficients. Solution 
of f and g can be expressed as 
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f =- !i_ s= w 1 
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* i/l+e ^ 


(2sin3/k - sin^B) „ . 


TG S 


w 


v/l+e^^ 

w 


( 21 ) 


and 


u 

g = - = 


TG S« e S- 
2 w 3 


^2f2 


l_ (2sinB/k - sin B) g 


2fl°r' 


( 22 ) 


VThere S^, S^, S^, and are convergent series given by 


a> 4p-3 °° 4p-2 oo 4p-l °® 4p 

c = y H q = y ■ . , ■ q = y ^ q — y IL— 

1 (4p-3): * ^2 (4p-2): » ^3 (4 p-d: ’ ^4 4p: 

P=1 p-1 P“1 p-± 


(23) 


and ri = y;o<r)<l. riisa dimensionless variable which gives 

the maximum value of distance from the wall (y"^ ) up to which the 

max 

solution is applicable for a given value of the angular velocity of 

rotation. For a typical flow in turbomachinery, 10 r' 10^. 

The approximate orders of magnitude of the first four terms are 1, 

+ -4 -7 

5x10 y/e,3xl0 y /e,8xl0 TGy /e for radial velocity 
w ' w „ w -5 

Q ■ > —7 

profile and 1,5x10 TGy,3xl0 ey ,8x10 y for stream- 

w 

wise velocity profile respectively. 

Very near the wall, it is sufficient to retain only the terms of 
2 3 

the order r)> 0 » and 0 . Hence the velocity profiles near the wall can 
be represented by 




+ 0 
w^ _ (2sinB/k - sin B) y , 

T ■ " 

l+e 
w 






(24) 


l+e 
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( 25 ) 



Thus, the velocity profiles in the inner layer are functions of 3, TG, 
k, Sf* and r* . 

The results from Equations (24) and (25), for k * 1.0, 3 = 0°, 30°, 
45°, 90°, J2° = 10 r* = 10^^, TG = 0, ± 1.0 plotted in Figures 4 and 
5 indicate that the flow is collateral, approximately up to y"^ = 10 or 
u/u^ = 8. When 3 = 0 (i.e. axial blading), both the radial and stream- 
wise velocity profiles are collateral for zero chordwise pressure 
gradient (TG = 0.0) and the velocity profiles follow the three- 
dimensional solution in stationary boundary layers with zero streamwise 
pressure gradients. For 3 > 0°, the flow is not collateral. The 
departure is large for radial velocity profile and increases with 3» 
reaching a maximum for 3 == 90° (flat sector) . Effect of chordwise 
pressure gradient parameter (TG) is large on the streamwise velocity 
profile and negligible on the radial velocity profile. The radial 
pressure gradient parameter (k) has the opposite effect. The departure 
for the radial velocity profile is large and increases with a decrease 
in k and is negligible for a streamwise velocity profile. 

The effect of rotation parameter Q° is to increase the radial 
velocity and decrease the streamwise velocity, compared to the no 
rotation case. It also has a large effect on the thickness of the 
inner layer. For lower values of U°, the inner layer is thicker. The 
effect of r* , which varies directly with Reynolds number is to reduce 
the departure from collateral behavior. In the limit r* Increases to 
infinity, the flow becomes collateral. 
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Figure 4. Velocity Profiles in the Laminar Sublayer of the Turbulent 
Blade Boundary Layer and Laminar Blade Boundary Layer. 



Figure 5. Hodograph Plot of Velocity Profiles in the Laminar Sublayer 
of the Turbulent Blade Boundary Layer and in Laminar 
Blade Boundary Layer. 





It Is thus evident that the collateral region exists only very 
near the wall for flows with radial or chordwlse pressure gradient. 

This Is true for the boiindary layer on both rotating and stationary 
channel. In the absence of these pressure gradients, the collateral 
regions extend farther from the wall. These conclusions, unfortunately, 
cannot be verified experimentally since measurements very close to the 
wall are difficult to make and are Inaccurate. 


2.3.2 Outer Layer (Region III - Figure 2) 

+ 3 

In the outer part of the boundary layer (l.e. for y > 10 ) , both 
the turbulent and viscous stresses In Equations (7) and (8) are small. 
Correlation between the two velocity components Is mainly governed by 
Inertia, pressure and Coriolis force terms In the momentum equations. 
Neglecting the shear stresses, eliminating v from the two equations and 
writing = f(u/U^> (26) 

the resulting equation can be written as (the procedure Is similar to 
the one carried out by Lakshmlnarayana et al. (17) ) 

f£’ {g- (sin^B+ RG) - £.{1 _ + f TG 

e e e 


- f^RG + stnh (1 + - 1 - - 0 (27) 

e e 

R ^^e 

where RG = ~ Is a radial pressure gradient parameter and R - r/r 
is the non dimensional radius. Equation (27) Indicates that. In general. 


f Is not only a function of ~ but pressure gradient parameters RG and 
TG also. For = flr (or k ■ 1.0) and sln3 ■ 1» the ^ove equation Is 
satisfied, if 



( 28 ) 
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even in the developing region of the turbulent boundary layer with free 
stream pressure gradients in x-*direction. This expression is identical 
to that derived by Senoo and Nishi (26) for a three-dimensional 
turbulent boundary layer in a stationary curved duct. 

I^en TG » 0 (i.e. no chordwise pressure gradient) which is nearly 
so for a fully developed channel flow (in the absence of interference 
effect near the outer and inner walls) or for very lightly loaded 
blading. Equation (27) can be simplified to 


g{2RG/(RGfsln^6)} , U ^ ^kslng " ^ ” U ^ 

gt + ■■■ y - • — (29) 

^ 2 slnB _ u j (RGfsin $) 2 sinB u j 

k(RGfsin^3) k(ROfsin^$) ^e 

2 

Where g = f . Since k, RG and 3 are now functions of r only, the 
above equation can be solved to give 


g = = C^[(^) 

e 


2 sin 3 j (2RG/RGfsin^3) 
k(RG+sin^3) 


+ (|” - sin3) 


sin3 _ y U X 2 
RG ^ 


(30) 


e 

where is a constant. The value of constant is obtained from the 
fact that solution (30) must be identical to general solution (28) when 
sin3 =k=li.e. C^=0. Hence the cross flow profile in the outer 
layer of the fully developed boundary layer in the channel is given by 




sin 3) 


sin3 /JL\ ^ 
RG ” 

e 


(31) 


Thus the cross-flow in the outer layer is not a simple function of the 
streamwise velocity, but of pressure gradient parameters k and RG and 
the blade stagger angle 3 also. It should be pointed out that k, RG 
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and 3 Inter-related by the local velocity triangle for a machine 

(as shown in Figure 1) and therefore can result In radial outward or 

Inward flow according as the centrifugal force Is larger or smaller 

2 

than radial pressure gradient (- Vq /r) respectively. 

In cases where the radial flow changes sign In a boundary layer, 
(e.g. k > 1 or when (2/k - sin3) sin 3/RG < (u/U^)^), Equation (30) is 
not valid, because the turbulent shear stress gradients are not 
negligible when the cross flow velocity changes Its sign inside the 
boundary layer. This has been confirmed by turbulent shear stress 
measurements in the Interference region reported in Chapter V where the 
radial flow velocity changes sign (i.e. from outward to Inward) and 
turbulent intensity and stress gradients are not negligible. 

Equation (31) indicates that the radial flow in rotating turbulent 
blade boundary layers increases if the factor (2/k - sin3) 

KG 

increases and is zero only in case k = sin3 = 1.0 I.e. when the centri- 

2 

fugal force is exactly balanced by the radial pressure gradient (Vq /r) 
at the edge of the boundary layer. 

2.3.3 Interference Layer 

Near the blade tip the rotation induced outward radial flow Inside 
the blade boundary layer is deflected by the annulus wall, producing 
an inward radial flow. This region is designated as the "interference 
region" and is shown in Figure 2. The extent- of this region Increases 
downstream of the blade passage, particularly when the flow becomes 
fully developed, resulting in a large decrease of streamwlse velocity 
component u. The flow in this region is complex due to the interference 
effects of the annulus wall and adjoining blade boundary layers. 
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The Intense mixing resulting from the Interaction produces large shear 
stresses and so viscous and turbulence effects cannot be neglected. 
Turbulent stress measurements In this region reported In Chapter V show 
that the stresses are not constant but vary linearly with distance from 
the blade surfaces. Based on these observations, an asyn^totlc solution 
Is derived below. 

Since the flow Is assumed to be fully developed, the velocity 
components u, v, and w do not vaiy with x, l.e. 

u = u (r, y) 

V «» v (r, y) (32) 

w = w (r, y) 

For the case of an Inducer, 3 - 80® - 84° (l.e. sin3 - 1). The equations 
of motion (7) to (9) with the uisual boundary layer assunqptlons and for a 
fully developed flow are simplified to. 


, 3u , u , 3u 1 X 1 ^ 

3r r 3y p dy p 3x 


(33) 


^ 3w 1 v2 1 ^”^r 1 ^ 

w T — I- v r (i2r - u) = - T a! T 

3r 3y r p 3y p 3r 


(34) 


3r r 3y 


(35) 


The left hand side of Equations (33) and (34) are functions of r 
and y only; therefore, the right hand side must also be Independent of 
X. Using the assunq>tlon that turbulent shear stresses are linear with 
y, and pressure variation with y in the boundary layer can be neglected. 
Equations (33) and (34) can be simplified by differentiating partially 
^ith respect to y and using the continuity Equation (35) to give, 
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3r9y r 3y 


3y 


v2-a. 0 

3y^ 


( 37 ) 


vhere F - flr - u. 

Equations (35), (36) and (37) can be solved by the method of separation 
of variables by writing, 

w = A(r) B(y) 

Fr C(r) D(y) (38) 

V - E(r) F(y) 

The resulting equations are: 

0 (39) 

0 (40) 


AB {C*D* - — C D'} + A B'C’D - A*B C D* + E F C D" 
r 


AB {A*B* - “ A B'} - 2 ^ C D* + E F A B" 


. . A* + A/r. - 

i- ( g > ^=4 


(41) 


Here ( ' ) refers to total derivative of a function with respect to its 
variable. The Interference Flow Regime for the Inducer channel Is 
shown In Figure 6. The tip clearance and leakage flows are neglected 
and the annulus wall Is assumed to rotate with the blade. The action 
of the annulus wall Is assumed to deflect the outward radial flow Jet 
near the blades, by 180° Invlscldly, to the main stream at the center 
of the section. The Jet mixes with the flow and penetrates up to the 
radius r^. The boundary conditions are: 


8w 8u _ n. 

^ - V = 0 at y = s 


(42) 


w “ 0 


at r “ r 
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Figure 6. Cross-Section of Inducer Channel Showing 
"Interference Region". 


3 . 


2S is the channel width as shown in Figure 6, The boundary conditions 
(42) in terms of new variables (Equation (38)) are: 

B* = D* = F = 0 at y = s 

and (43) 

A = 0 at r = r 

c 

Equations (39) and (40) can be separated only if 

B D* = c.B* D = c, f D" (44) 

D O 

and 

B B' = c_D D’ = CoF B" 

/ o 

where c_, c,, c_, and c^ are constants. 

5 6 7 8 

Equation (44) shows that 

D - Cg B (45) 

where c^ is a constant. Also 

2 

c^ = 1, = Cg, c^ = 1/cg (46) 

From Equations (41) and (44) , it is clear that, 

CqC, B B" + B B’ Itt = (B B" + B*^) (47) 

O 4 D 


It is assumed that B can be written as a polynomial series of y as 
follows : 

B = b^ + b^ (y-s) + b^ (y-s)^ + bg (y-s) ^ + . . . + b_^ (y-s) ^ + 

(48) 


where b , b_ , b_, b _ are constant coefficients. Since B is finite at 
o 1 2 — 1 

y = s, b_ 2 » . . . b_^ terms are all zero. Also B is symmetric 

about y = s, hence 


b 


3 



“ ^2n+l 


= 0 


( 49 ) 


where n is an integer. 


Using the boundary condition B'(y=s)=0, 


( 50 ) 


Therefore B can be reduced to the following sliq>llfled form: 


B ■ % + ^2 *4 • * • + + . • . 


2n 

where x ■ • Substituting for B from Equation (50) in Equation (44), 

following relation between coefficients is obtained: 

2 ,2 ^ 

+ X + . . . H 

m=0 


+ . . . + Z (2n+2-2n) 


c.c„ [+2a„a,+(-r^+12a„a,)x^+ ... I — 


o 4' 
2 


'4 8 


o 2 ' 3 


o 4 


m=0 


where m - 0, 1, 2 . . . etc. 

Comparing the coefficients of x of same index. 


^2m^2iH-2-2m^ 

(51) 


2n 


^^4^8 “ ^ 


2a, 


2 ““2 
2a« + 4a a, = ■ ■ ■ - + 12 a a, 

2 o 4 3 o 4 


For general case 


Jo (2»f2-2») a2^a2^2_2^ =J^ ~ ^2m^2i^2-2m 


(52) 


(53) 


This is simplified to give 


n 


“ (2n+2-2m) (2n~4m) ^ . 

(2mfl) ®2m^2n+2-2m 

Using Equation (54) all the coefficients of B can be found, 


(54) 


2 3 4 

1 /^2 ^ 1 /^2 s 1 /^2 

®4 “ ? ^a ^ " 90 ^ 2^ ’ ^8 “ 2520 ^ 3^ * * * 

o a a 

o o 

The solution of B in Equation (50) can be written as 

*2 2 ^2 2 ^2 3 Y^ 

B “ X' + (/)" (/> lo + • • •’ 

O O D 


(55) 
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Substituting for B In Equation (45) , D becomes 

^2 2 ^2 2 

D = c. B = c^a [l+(^) X + (^) I-+ . . . ] 


'9 o" ^a 


a ' 6 

o 


( 56 ) 


Solution for F Is obtained by Integrating Equation (41) with respect to 
X« using the Equation (55) for B and the boundary condition F « 0 at 
X “ 0. 

F = c,a [x+ (^) (^)^ + . . .] (57) 


'4 o' 


‘a ' 3 'a ' 30 'a ' 630 

o o o 


Substituting for E from Equation (41) and using Equation (44), the two 
Equations (39) and (40) are reduced to 


and 


AP 

AC* - ~ + AC* - A*C + (A* + A/r) C =• 0 


AA* - “ ‘ \ + (A* + A/r) A “ 0 


Equation (58) gives 2AC* =0. Since A 0, 

■ <=10 

where c^^ is a constant. Equation (59) gives, 
2AA' -(2/r^)C^ c ^ = 0 


(58) 


(59) 


(60) 


(61) 


Substituting for C from Equation (60) and integrating A, is obtained as 

(62) 


a2 2,2^ 

^ ^ - ^10 "^9 


where c^^ is an Integration constant. Since A =• 0 at r =» r^, A is 
given by, 


A = c« c 


/ (-^)-(4r: 


9 10 ^ 2 ' ^ 2 

r r 
c 


(63) 


Substituting for A from Equation (63) in Equation (41), E is obtained as 
follows 

E = {(r/r )^ - (64) 

r ^ 

c 
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The solution of u, v, and w Is obtained by substituting the expressions 
for A, B, C, D, E, and F from Equations (63), (55), (60), (56), (64), 
and (57) respectively as follows: 


W_ . 1 r ^^~^cV /2 a+A, 2 .^2, 2 f2,3 , 

w R \ „ 2 ^ ^ X + la J 6 ^ ^ 90 ^ 

t 1-R o o o 

c 


(65) 


t t o o o 


( 66 ) 


y-= )sl+ (^)2 £+ (f2)3 xL+ . } (67) 

o 1-R o o o 


where w^ and u^ are the radial and chordwlse velocities at the tip at 

2 2 

y * s (l.e. X = 0) • ^ constant. The ratio (a^/a^) = (3 w/3x /2w) 

at X = 0 ^ negative constant because w Is negative and has a minima 

at y = s. Its value Is to be determined from the experimental data. 
From Equations (65), (66), and (67) the following correlations between 
u, V, and w are obtained; 

w 


w 


f2r-u 


2 2 
^ R -R ^ 
t N r C 1 1/ 2 

O J 


= ( — ^){- 

Qr -u 1-R ' 
t t c 


R^-R ^ , /2 
= e { 

max 1^2 


1-R 


( 68 ) 


and 


V R 1-R 0^0 o 

2. = _2_ { £_}(2::S.){1 _ 2 (_2. ,2^,2 y 

w w. „2 „ 2^^ s 3 2' ^ ' 

t R -R os 

c 


(69) 


where e Is the maximum radial Inward flow angle at the tip. The 
max 

results of Equations (65-69) plotted In Figure 7 show the following 
properties of the velocity profiles In the Interference layer: 
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(Eqn . 66) 
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2 — (Eqn. 65) 
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Velocity Profiles at R = 0.98, x/c =* 0.9, 
Open Throttle Data in Figures 26 and 32) . 
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Hodograph Plot of Radial Velocity Versus Streamwlse Velocity 
Defect. 


Figure 7. Velocity Profiles in the Interference Region. 



1. The radial Inward flow is maximum near the tip and decrease to 

zero at a radius r^. Its profile Is an even polynomial in terms of 

distance (y-s) from the center of the channel and Is approximately 

2 2 

parabolic. The ratio (* 2 / 3 ^) Is given by (3 w/3x /w) at X ■ 0 (l.e. 
y *■ s) and Is directly related to ratio of shear stress gradient to the 
magnitude of velocity In the radial direction at the center of the 
channel. 

2. The velocity defect in streamwlse velocity varies directly 
with radial velocity at any radius. 

3. The velocity normal to the blade v Is zero at mid channel, 
Increases with distance away from centerline and Is directly proportional 
to the product of radial velocity and distance from the center of the 
channel. 

2.4 Momentum Integral Analysis 

In the field of turbomachinery. It Is In^ortant to know the gross 
boundary layer growth characteristics In a rotating passage. In 
particular, for a turbomachinery designer It gives an estimate of 
aerodynamic losses, regions of separation (If any) and the severity of 
the three dimensionality of the flow. The earlier analysis In section 
2.3 Is aimed towards an understanding of the structure of the boundary 
layer. This has both scientific and engineering relevance, whereas the 
Immediate need of a turbomachinery engineer Is to know the boundary layer 
growth parameters viz. the momentum thickness, the limiting streamline 
angle and the shape factor, which are necessary In predicting the 
overall aerodynamic performance of a turbomachine. The following 
momentum Integral analysis enables the calculation of turbulent boundary 
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layer growth in a rotating blade passage. The chief virtue of an 
integral method for turbulent boundary layers lies in the implicit and 
global manner in which the effects of turbulence are Incorporated. 

The method la an extension of Lakshmlnarayana et al.*s (17) method 
for predicting the boundary layer growth on a single* rotating helical 
blade. It includes the effect of the pressure gradients iijq>osed by the 
Invlscid external flow. The two momentum Integral equations (one in the 
mainstream and the other iii radial flow direction) are reduced to two 
first order partial differential equations using a power law profile 
for the mainstream velocity, Mager's profile for the i cross-flow and a 
skin friction relation based on the experimental data of Anand and 
Lakshmlnarayana (27) , The skin friction coefficient is a modification 
of Ludwig and Tlllmann*s (28) skin friction correlation for two- 
dimensional flows to three-dimensional turbulent boundary layers with 
pressure gradients. Including the effect of rotation. The exponent of 
the power law for the streamwise velocity conq>onent is a function of the 
local shape parameter (H) , which is defined as the ratio of streamwise 
displacement thickness to streamwise momentum thickness. It thus 
represents the shape of the streamwise velocity profile. 

In turbomachinery blade passages, the flow experiences large pres^ 
sure gradients. These pressure gradients cause large changes in 
velocity profiles and consequently in the shape parameter H, The 
variation of H cannot be neglected and an additional equation is required 

Out of the available auxiliary equations described by Kline et al. 
(29) , only the "energy Integral equation" and the "entrainment equation" 
have been found suitable for the turbulent boundary layers. The former 
involves more empiricism in calculating the dissipation integral term. 
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and assumes the shear stress distribution in the turbulent boundary 
layer. Since the shear stress variation differs largely from flow to 
flow, especially those under the influence of rotation and curvature, 
the entrainment equation rather than energy integral equation is used. 
The entrainment equation is derived from the concept that a turbulent 
boundaiy layer grows by a process of "entrainment" of the inviscid fluid 
at the edge of the boundary layer into the turbulent region and has been 
found to give good predictions in a variety of flows as demonstrated by 
Kline et al. (29) . The momentum equations and the entrainment equation 
are solved numerically by a finite scheme to predict the three most 

important boundary layer parameters viz. 9 , t > £ > and H, of the flow. 

11 w 

The following two momentum equations in rotating coordinate system 
(x, y, r) can be derived from Equations (7) to (9) in a manner similar 
to Mager (30) , 

99. . . 9, 


^®11 ®12 
+ (2+H) TG + (2 
dx r r 


*^ 2 * “'^12 2 '^12 
“^) RG + + (1+sin 3) 

r 9r ' r 


^2* 2 ^fx 

+ sin 3 (^ - sin3) = 


(70) 


3(r922) 

r3r 


9 


+ 2 


22 

r 


®21 ^®21 ® 11 ^^^ ^ 

RG + 2 — TG + -- — - + 

r 9x r 


sin g _ sin3) + £ ~ 0 (71) 

r k w 2 

where 9--, etc. are momentum thicknesses, <5 *, d * are displacement 

thicknesses in mainstream and cross flow directions defined as follows. 
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11 


c 

■I 


e 


^)dy, 

e 


(1 - 


^)dy 

e 


'12 


52 - (1 

U 

e 


^)dy, 

e 


c 

■f 


2- dy 
U ^ 
e 


(72) 


'21 


o 

I 

“ u2 


dy 


022 = 1 (^)'dy 


where C- = 2 t /pU , T is the streamwlse component of surface shear 

X3C 3C 0 3C 

o o 

stress, T = e T (the flow is assumed to be collateral near the 
r w X 
o o 

surface) and k = U /Qr, U is the velocity in mainstream direction at 
e e 

the edge of the boundary layer and = skin friction coefficient. 

The assumed mainstream velocity profile is 




(73) 


e 

where 6 is the local boundary layer thickness, defined as the normal 
distance from the blade surface, where the strearawise velocity reaches 
its first maxima. H is the local shape factor and n is the 
non dimensional y coordinate. The assumed cross-flow profile is (Mager 


(30)) 

|-e^(l-r,)2 (74) 

where is the limiting streamline parameter (tan o^) as defined in 
Figure 3. 

Equations (73) and (74) simplify Equations (70) and (71) in terms 
of four parameters H, and The resulting equations are: 
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ae.. 9-- £0 ry a 

+ (2+H) TG + J (1+sln 3+2RG) + ^ 

+ - - sin^3- RG) — ~ = -|^ (75) 

and 

2 

“tkr ^ ^ •" ®ii ^ ^ -4^ 

+ ^7 w®ll«> > - H(^ - sln^B) ^ - 0 (76) 


where 

* 

_ ^12 (30+14H) ^ _^2 16H 

^ ^w®ll * ^w^ll ” 


(77) 


®22 24 „ ^21 2 

^ " £ 0,/ “ (H-1) (H+2) (H+3) (H+4) ’ ^ " £ 0^ , “ (H-1) (H+2) 
w 11 w i± 

The skin friction coefficient C_ is related to the other three 

fx 

parameters 0^ - , £ » and H on the basis of experimental data by 
11 w 

Lakshminarayana et al. (17) and Anand and Lakshminarayana (27). In 
flows with zero or mild pressure gradients the skin friction correlation 
based on the experimental data of turbulent boundary layer on a single 
rotating blade by Lakshminarayana et al. (17) is 

= 0.0582 R^ (1 + 0.85 /(£^x/r)) in the developing (78a) 
region and 

-1/5 

= 0.079 R^ for fully developed region (78b) 

where R^ and R^ are Reynolds numbers, both based on the local streamwise 
velocity at the edge of the boundary layer and distance x from the 
leading edge for R^ and radius r for R^ respectively. The skin friction 
relation for flows with large pressure gradients is based on the 
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experimental data for a turbulent boundary layer flow In a rotating 
channel, described by Anand and Lakshmlnarayana (27) , 

‘=fx ■ “l >*8^^ “ (1+B^^e^(x-X^)/C ) (79) 

where a_, b, Y_ and 3^ are universal constants. R-, is the Reynolds 
111 ®ii 

number based on the streamwlse velocity at the edge of the boundary 
layer and streamwlse momentum thickness is the chordwlse 

distance of the transition line from the leading edge. 


2.4.1 Two Equation Model - Mild Pressure Gradients 

For zero or mild pressure gradients in the streamwise direction, 

the shape factor H, may be assumed constant (H = 1.286) and the number 

of variables are reduced by one. The two momentum equations together 

with the assumed skin friction relations 78(a) and 78(b) are sufficient 

to solve for 0^ , and e and there is no need for an additional equation 
11 w 

for H. Equations (75) and (76) using from Equations (78a) and (78b) 
and Equation (77) with H = 1.286 are reduced to; 

~+ 0.542 e |^+ A {3.285 — + 0.543 + 

dx w dr r <3r 


_w „ 1,59 RG - 2.134 sin^3) + 0. 335-0. 2/x}-.0296A « 0 

r k 


and 


( 80 ) 


e ■—+ 0.530 £ ^ |^+A{-^ (0.318 + 1.05RG)+ e (^ - — 
w 3x w 3r r ' w r x 


0£ 8e j 

+1.059 TT-^)+ (0.715 sin 3“ 0.8sin3/k)} + 0.014e A=0 

dr dx w 


(81) 


where J » 0.542, L = 2.673, M = 1.129 and N * 2.13 are used in the above 
equations, for H ^ 1.286. 


45 



For the present case of inducer, the radial coordinate is 
normalized by tip radius (r^ = 1.5 feet). The above equations are 
further simplified to, 

3A . . TG 


t ^ 0.362 t -^A{2.19 0.362 + Ji (1:^ 


sin3 


and 


-1.06RG-1.423sin 6) + 0. 223-0. 2/x}-0.0296A * 0 


e |^+ 0.353 e ^ A{ (0 . 212+0 . 7RG) -^ + e (1.33 

w dx w dR R w R 


( 82 ) 


0.2 

X 


+ 0.706 -^)+ '^+(0.715sln^3 - ^ - ) }+0 . 014e^ A = 0 (83) 


where 


1/5 

A = 0., , R is the momentum thickness parameter, and A = 1+0.85 
11 X 


/(£ *x'/r) is a rotation parameter which gives the direct effect of 
w 

rotation on skin friction as explained by Lakshminarayana et al. (17). 
The boundary conditions are: 

1. At the leading edge, A(0,R) = e (0,R) = 0 

w 

2. At the hub, A(x,0.5) and £ (x,0.5) are determined from single 

w 

blade solution which is modified to include pressure gradient parameters 
RG and TG. 

The flow near the leading edge is found to be laminar and its 
extent (x^) has been determined by experiment. Banks and Gadd’s (7) 
analysis is used to predict 0^^ and in the laminar region up to the 
point of transition. It is to be noted that this analysis is valid 
for the laminar boundary layer over a single rotating blade, and so does 
not take into account the effect of imposed pressure gradients on the 
boundary layer. Since the actual pressure gradients in this region 
are very small for the present case, neglecting this effect may be 
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justified. For x > x^, the flow is assumed to be turbulent and the 

turbulent flow Equations (82) and (83) -are used to predict 0,- and e , 

11 w 

The partial differential Equations (82) and (83) are solved 
numerically using fourth order finite difference scheme. The values of 
TG, RG are assumed to be known from the inviscid analysis of the channel 
flow. For the case of rotating helical channel used in this investi- 
gation, these values are derived from the measured blade pressure 
distribution. The convergency criteria used in the numerical analysis 
are 

< lO"^ feet 

(84) 

- e5”7^^| < 10“^ radians 
i» J i » j 

where n and n-1 are nth iteration and (n-1) iteration respectively. A 
zero order solution is obtained by solving for A from Equation (82) by 
neglecting terms containing and its derivatives. This is used as 
initial data for Equation (83) to solve for and thereafter Equations 
(82) and (83) are successively iterated till convergency, according to 
the criteria of Equation (84), 

The numerical solution of flow over a single rotating helical blade 

I 

t 

was also obtained and found to agree with the analytical solution 

obtained by Lakshminarayana et al. (17), where radial variation of e 

w 

is neglected. The radial variation of was found to be negligibly 
small for a single blade and justifies its neglect for that case. 

2.4.2 Three Equation Model - Large Pressure Gradients 
For botindary layers with large irqposed pressure gradients which 
normally occur in a turbomachinery operating at the design and off 
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design flow conditions, variation of shape factor (H) is large and an. 
additional eqxaation to account for its variation is required. The 
entrainment equation, which as described earlier is better than other 
auxiliary equations (Kline et al. (29)) and has the least empiricism, 
is derived below. 

The entrainment equation is obtained from the integration of 
continuity equation from blade surface (y = 0) to the edge of the 
boundary layer (y = 6) . The continuity Equation (9) is the same both 
in the rotating and stationary coordinate system (x, y, r) . 

Following Nash and Patel (31), its integral form can be written as 


h ^ ^ If ^ «e f - 


(85) 


where 6^ and 6^ are displacement thickness parameters in strearawise 
and cross flow directions respectively and are defined in Equation (72) . 
E on the right hand side represents the volume flow rate per unit area 
through the surface 6(x, r) and is the rate of entrainment of inviscid 
external flow into the boundary layer. 

The entrainment process is a highly complex and inherently unsteady 
phenomena and its direct measurement is difficult. The entrainment rate 
depends both on the mean flow parameters, such as the streamwise 
velocity defect, the rate of boundary layer growth, streanwise velocity, 
displacement thickness etc., as well as on the turbulent quantities such 
as turbulence intensities and stresses in the outer layer. The stream- 
wise velocity defect in turn is directly related to streamwise velocity, 
the shape factor, the magnitude of cross flow and the blockage effects 
due to confinement of external inviscid flow in the channel. 


48 



The effect of streaiwlse velocity and the shape factor on the 
entrainment rate have been extensively studied, resulting in Head's 

(32) entrainment equation, 

E = *> (86) 

where ^ * is a shape parameter, defined by ^ * = (6-6^^ ^^^11 

uniquely related to the shape factor H. For one parameter main velocity 

profile (according to Head (32)) the equation for shape parameter is, 

-2 71 *5 

Hg^g * - 1.535(H-0.7) 3.3 (87) 

and 

*) - Cg (H^_^ * - (88) 

where C_ is the entrainment coefficient. Its value for two-dimensional 
E 

boundary layers is 0.0306, Qualitatively, it is known (Kline et al. 

(29)) that the effect of an increase in turbulence level in the external 
stream is to increase the entrainment rate. On the other hand, an 
Increase in cross flow causes a reduction in boundary layer thickness 
and hence its growth rate in streamwise direction. Since the entrainment 
rate is directly proportional to boundary layer growth rate, consequently 
an increase in cross flow causes a decrease in it. Assuming that the 
variation of entrainment rate with H^_^ * follows the same relationship 
even for three-dimensional flows, e. g. Cham and Head (10), Lewkowicz 

(33) , the effects of turbulence level and cross flow can be incorporated 
by assuming that C_ = C_ (q , e ) , where q is the free stream turbulence 
level. Hence the total entrainment rate can be written as follows: 

E - U^Cg(q^. - 3.0)-°- 658 (89)- 
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Substituting for E in Equation (85) , the auxiliary or entrainment 
equation becomes 


-0.653 


{U^(6-6^ )} + ^ {U^r 62 }= 

Writing in terms of three variables 0_-, e , and H, the equation becomes 

li w 


|-{U . *)+ -|- {U Lr0 

9x e 11 6-6^ r dr e 




The entrainment Equation (91) is further simplified to give a direct 
equation for the shape parameter ^ , 


■ L. + u , *[^+ ^ (U ^ *-3.0) 


9x 


^6-6i ^r ’ 6^^ 9x ^ 6^^ '^"‘^1 


1 r 

7 [e„(l+RG + -5- 


X . 3ei 
3r 37^ 


(92) 


The two momentum integral Equations (75) and (76), together with the 
expression (79) for skin friction and the entrainment Equation (91), 
constitute a set of four equations for four variables ^ ^ 

can be eliminated from (75) and (76) by substituting (79) and these 
equations are simplified to: 

|^+ J£ |^+ - [TG{(2+H)(l+b)-b} + (l+b)r + 

dx w dr r or 


e { 2L(l+>) . sln6 _ RG(L(l+b)-J(2+b)-sin^e(L-J)(l+b)+J(l+b)}] 
w k 

“l * 

- A (1+b) = 0 


( 93 ) 



and 


e {(l+b)+RG(2+b) } +e {TG(2+b) + 

wdxNwdr r wN w 


2 dn\ (1+b) (1+H) , . 2„ 

N " 5 T~^ + “n + ~N ® 


a (1+b) * 

+ e A > 0 

2N w 


2H(l+b) sin6 , 
N k ^ 


(94) 


where A = 6- , Rq is a new variable for momentum thickness and 

11 


A = 10 ^ (1+3t /e /(x-x )/e ) 

1 W _ t ' 

The ex^K^r-'' mental data of Anand and Lakshminarayana (27) indicate that 

the values of b = +0.268, = 0.678, 3^ = 0.52, = 0.172. Equations 

(93) and (94) are to be solved siraulatneously for 9^^ and e^, using the 

same finite difference scheme as described in the section 2.4.1 and 

assuming initially a constant value of H. These values are to be used 

in Equation (92) to find a new value of H. Thus Equation (92) is to be 

solved successively with Equations (93) and (94) till the convergency 

criterion (84) for A and together with the following criterion for H 

w 

is reached. 

I ^ 0.001 (95) 


where and are values of the shape factor H at nth and (n-1) th 

iterations respectively. The boundary conditions and the procedure for 
incorporating laminar solution are same as described in section 2.4.1 
Since an expression for C for three-dimensional flow is not yet 

£r 

available, no numerical solutions are given in the thesis, even though 
the computer programming has been conqileted. 
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III. TEST FACILITY AND DATA ACQUISITION EQUIPMENT 


3.1 Experimental Test Facility 

All the measurements reported In later chapters were performed 
using the NASA Axial-Flow Inducer Test Facility, shown in Figure 8 and 
located in the Turbomachinery Laboratory of the Department of Aerospace 
Engineering at The Pennsylvania Sta^e, University . A 5 hp Poly dyne 

' •* ■ • • • • '-' 4 ' > 

Adjustable Speed tfotor (speed ranged 286-2000 rpm) drives the rotor 
through a belt and pulley system. M aero dynamic ally designed throttle 

* . ■ - ' i' ' ' ' ' , 

can vary the flow rate through t|ie machine, so that the rotor can be 
Operated at different flow coefficients on the performance map. The 
annulus has a well designed bell mobth inlet so as to give uniform 
axial flow at the inlet to the rotor. The facility has a number of test 
windows and static pressure holes so that a comprehensive flow survey is 
possible. Extensive flow visualization and quantitative measurements of 
the flow field were performed upstream, inside and at the exit of an 
axial flow inducer rotor (referenced 17, 21, 22, 23, 34) by varying the 
number of blades from 1 to 4 and at different flow coefficients. These 
measurements would be compared to the data obtained with the four bladed 
flat plate inducer described below. 

Rotor Facility . The rotor consists of four equally spaced 
transparent plexiglass helical blades of constant thickness (0.5 inch), 
an axial advance of 10 inches and an angular extent of 300 degrees. The 
blade stagger angle varies from 78^ at the hub to 84^ at the tip. Other 
dimensions of the rotor and facility are shown in Figure 9. Both the 
leading and trailing edges of the blade are shaped to form a symmetric 
rfoll shape to ensure smooth flow at both the edges. The blade 


52 




0 ANNULUS WALL 


i.r 

t34«iaii 


10 “ 

(25.4 cm) 



Figure 9. 



surface can be considered hydraulically smooth in the operating range 
because <2.5 where X is the mean height of the surface roughness. 

To facilitate blade static pressure and limiting streamline angle 
measurements, the adjoining blades of the channel have ten static 
pressure taps at each of the five axial locations. The rotor is 
statically and dynamically balanced before installation in the Inducer 
test facility to reduce vibrations. 

One unique feature of this facility is the traverse unit for the 
rotating probe measurements. Referring to Figure 9, it consists of a 
light traverse disk driven by a fractional horse power step motor. The 
probe holder and the probe is mounted on the traverse disk. The probe, 
which is located inside the passage, traverses in a slot in the hub. 

The traverse disk and the probe can be traversed circumferentially in 
steps of 1.8° by the stepping motor, thus enabling accurate location and 
measurement of the flow characteristics inside the blade boundary layers. 
The motor is powered and controlled by a stationary traverse indexing 
device through an eight channel slip ring unit. The traverse unit and 
the probe are locked* except when indexed. Thus the circumferential 
traverse can be accomplished, through the remotely controlled index 
device, without stopping the rotor. The radial movement of the probe 
is carried out manually. Once all the radial and circumferential 
traverses have been completed, the entire traverse assembly is moved to 
a different axial location and the measurements repeated. The instrument 
package, consisting of commerlcal slip ring unit, mercury slip ring 
unit, and a pressure transfer device is housed in a stationary cowling 
attached to the nose cone as shown in Figure 9. In order to eliminate 
possible transmission of the rotor assembly vibrations to the instrument 
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packages a flexible rubber coupling is used. All the measurements have 
been carried out at the rotational speed of 450 rpm corresponding to a 
Reynolds number of 6,7 x 10^, based on the tip radius. The speed of the 
Inducer is maintained constant up to an accuracy of 0,1 rpm by means of 
a photocell circuit with a 60 slot calibrated disk mounted on the rotor 
shaft and displayed on an electronic counter. For checking the probe 
position and for photographing the flow in rotor passages, when rotating, 
a high Intensity strobeflash unit is used. A schematic of the experi- 
mental setup is shown in Figure 10, The following devices are used for 
measurements in rotating system. 

3, 2 Data Transmission Devices 

Three Channel Pressure Transfer Device (ffTD^ . The pressure transfer 
device (PTD) has three channels, each of which is made airtight by using 
double sealed ball bearings. Pressure leakage along the cylinder is 
prevented by the use of 0-rlngs and the plastic sealers. Pressure from 
the rotating system is transferred to the pressure transverse device 
through a hollow shaft which is fastened to the nose cone, and is measured 
by a micro-manometer. In order to reduce the possible interference on 
the incoming flow, the pressure transfer device is housed Inside a 
streamlined cowling. 

Ammonia Transfer Device (ATP) . The direction of the flow near the 
blade surface is determined by means of the ammonia streak method, i.e., 
by releasing a small amount of ammonia gas through the holes drilled on 
the blade surfaces and measuring the resulting ammonia trace on an 
ozalid paper strip attached radially on the blade surface adjacent to 
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Figure 10. Experimental Set Up (Open Throttle Measurements) 












the measuring station. ATD transfers ammonia gas from the stationary 
system to the rotating system. The ATD is mounted on the stand outside 
the rotor assembly. 

Slip -Ring Unit . A commercial eight channel slip-ring unit is used 
to conduct electrical signals from the stationary reference frame to the 
rotating reference frame of the rotor traverse unit. Electrical 
continuity is provided by carbon brushes in contact with a rotating 
commutator aligned along the inducer's rotational axis. The slip-ring 
unit is mounted on the pressure transfer device, and all electrical and 
pressure connections were transferred through a hollow shaft and 
flexible couplings to the nose cone of the inducer. 

Mercury Slip-Ring Unit. A six channel mercury slip-ring unit 
(Rotocon, MSD /i^6, Meridian Lab.) is utilized in transmitting the hotwire 
signals from the rotating reference frame of the inducer to the 
stationary hotwire anemometers. The slip-ring unit exhibits the smallest 
and most stable resistance in the transfer of measuring signals from the 
rotating electrical elements to the stationary electrical conductors. 
Contact between the rotating wires and the stationary contact screws 
is made through a round contact disc to which the rotating wire is 
connected, rotating in mercury. Triple-distilled mercury is used to 
provide the greatest conductivity and the lowest noise level distortion 
possible. The slip-ring is mounted on a four arm vibration damper and 
thus eliminates any spurious electrical signals due to possible 
transmission of rotor assembly vibrations to its electrical contacts. 
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3,3 Instrumentation and Probes 


3.3.1 Measuring Equipment 

All the pressure sensitive signals are measured with a Merlam 
Micro -manometer. All the pressure measurements lie In the range 0.05 
- 1,1" of water. The minimum accuracy in pressure measurement Is 2%. 

Two dual channel constant tetq>erature hotwire anemometers are used 
to provide the three channel capability necessary for the three- 
dimensional turbulence measurements. The original anemometer circuitry 
has been designed by Lumley (35) to give a flat response (Appendix A) 
up to 50 KHz even for large turbulence signal levels (up to as high as 
20%) , Operating characteristics of the circuitry are described in 
Appendix A. The anemometer power supply is kept on for 24 hours before 
taking data to minimize the voltage drift. The number of ground loops 
with the connected circuitry are minimized to prevent surging in the 
anemometer circuitry and to reduce the electrical noise from power 
supply and other electrical equipment in the room. A schematic of the 
hotwire set up is shown in Figure 11. 

3.3.2 Probes 

A three hole probe disk type (Figure 12) is used to measure the 
total pressure and radial flow angle (defined as the angle which the 
local streamline makes with the streamline at the edge of the boundary 
layer) Inside the rotor blade passages in a measuring grid of 5 x 5 x 40 
in tangential, radial and blade to blade directions respectively. The 
probe is 1/4 inch in diameter and 1/16 inch thick. Three holes 1/32 
inch diameter each, drilled circumferentially on the disk and 45^ to 
each other form the sensing area. The disk probe has a number of 
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Figure 11, Schematic for Hotwire Measurements, 
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1/4” Diameter, 1/16" Thick Disk Sensors 
with 3 Holes Drilled 45° to Each Other 


1/16" 



Figure 12, Disk Type Botmdary Layer Probe. 
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advantages over a conventional yaw probe for measurements close to a 
blade. It can be easily aligned parallel to a blade surface^ produces 
less flow distortion and has a wide linear flow angle measuring range 
(±20^) as shown In the calibration curve (Figure 13). The calibration 
curve Is later corrected for total pressure and flow angle using the 
Open Jet Timnel Facility as described In Appendix B. 

Two Preston tubes of outside diameters 0.028 Inch and 0.065 Inch 
are used to measure the total wall shear stress by placing them flush 
with the blade In the direction of measured limiting streamline angle 
at that location. _ 

A special three sensor boundary layer hotwire probe (Thermosystems 
Inc.), shown In Figure 14, Is used In all rotating hotwire measurements 
to obtain the mean velocity components and the turbulent quantities from 

the direct measurement of the three mean sensor output voltages, 

— 2 2 2 

and the six products of fluctuating voltages » ®2 * * ®1®2* 

e^e^, ® 2 ® 3 * orientation of the sensors Is also shown In Figure 14 

with respect to the reference coordinate system x, y, r. Sensors 1 and 
2 are In the xy plane and are orthogonal to each other. Sensor 1 makes 
an angle with the x-axls. Sensor 3 Is located midway between 1 and 
2 and Is at angle bo the xy plane (l.e. the plane containing sensors 
1 and 2). The direction cosines, the dimensions of sensors and the 
support needles etc. are given In Appendix A. All six needles and their 
electrical connections are Insulated from each other so that each wire 
can be operated Independently of others. The probe body has an elliptic 
cross section (1/8" x 1/4") and therefore can be used for measurements 
as close as 1/16" from the blade surface without touching It. The experi- 
mental set up for hotwire measurements Is shown In Figure 15. 
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Figure 13. Calibration Curve for Disk Probe. 
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Figure 14. Three Sensor Hotwire Boundary Layer Probe. 





3,3.3 Signal Processing Equipment 


The output of the anemometer consists of three fluctuating a.c. and 

three d.c. signals. The d.c. signals are measured with a 3430A Hewlett 

Packard Digital Voltmeter up to an accuracy of 1 mv, through a three 

channel selector switch. The a.c. signals are passed through a 

multichannel selector switch. This switch has six input channels. Each 

two of these are connected to one of the three a.c. signals. The switch 

has two output channels and therefore gives any desired dual combination 

2 2 2 

of the three instantaneous signals (i.e. e^ , 62 » ® 1 ^ 3 * ® 2 ® 3 ^ ' 

The signals are now filtered through a dual channel 5 kHz low pass 
filter-amplifier circuit. Both the turbulent spectral measurements 
without the filter and the frequency corresponding to Kolmogorov scale 
(maximum of 3.5 kHz at x/c = 0.1) show that the signal is either of the 
order of or less than the electrical noise signal above 5 kHz and thus 
governed the choice of 5 kHz filter. The signals are multiplied in the 
Model 520 Analog Devices Multiplier Unit and integrated in an integrator 
up to an accuracy of 1 mv for a large integrating time (45 seconds for 
each observation) using a synchronized stop clock. The integrator is 
calibrated to have an integration factor of unity between input and 
output. The integrated output is read through a 3440 A Hewlett Packard 
Digital Voltmeter with a 3444 A Multifunction unit, up to an accuracy 
of 1 mv. A schematic of the hotwire set up for measurements of mean and 
turbulent quantities is shown in Figure 11. 

A type 55D35 DISA RMS unit with variable integration time is used 
for overall turbulence level measurements. A 302A Hewlett Packard Wave 
Analyzer is used for spectral measurements. The wave analyzer has a 
carrier frequency of 100 kHz and a bandwidth of 7 Hz for its narrow band 
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frequency selector filter. Its output is squared and amplified using a 
Weston Model 1056, Ballantine Laboratories Mean Square Meter, integrated 
for 60 seconds in the integrator and the final output is read on the 
3440 A Hewlett Packard Digital Voltmeter. Figure 11 also shows a 
schematic of experimental set up for spectral measurements , 

A four channel type 3A74 Four Trace Tektronix Storage Oscilloscope 
is used in all the hotwire measurements to visualize the magnitude of 
all three turbulent signals and their multiplied output. A Model 200CD 
Hewlett Packard Wide Band Signal Generator is used to calibrate the 
Wave Analyzer and to find the gain of individual components in the 
hotwire circuitry, except the hotwire anemometer channels, which are 
calibrated Individually in the calibration tunnel described along with 
other peripheral equipments in the following section. 

3.4 Peripheral Equipment 

A low turbulence calibration tunnel is used for the hotwire 
calibration. The horizontal wind tunnel has a square test section with 
each side of one and a half inch and operates within the range of air 
velocities of 0-300 feet per second. The calibration velocities are 
measured with a United Sensor 1/8 inch pitot-static tube and the micro- 
manometer described previously. 

Pressure sensitive probes (e.g. cobra probe, total pressure probe, 
etc.) are also calibrated in the same tunnel. 

A one-foot diameter open jet tunnel of the Garfield Thomas Water 
Tunnel Research Facility of The Pennsylvania State University is used 
to calibrate the three hole disk probe for total pressure and flow angle 
at different velocities. The calibration curves are given in Appendix B. 
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The local hotwire making facility enables sensor replacement In a 
comparatively short time. The probe Is mounted on a small table vice 
with sensor area facing an adjustable focus tilting microscope. The new 
wires are soldered using a special solder Iron and etched to required 
resistance by a dilute nitric acid solution. The desired resistance of 
the hotwire and the associated "over heat ratio of resistances" of the 
hotwire anemometer are measured up to an accuracy of 0.1 ohm using a 
null galvanometer. 

An optical profile projector of the Applied Research Laboratory 
Is used to measure the geometry of probe needles and the dimensions of 
hotwire sensors. It also measures the wire angles with respect to the 
reference coordinate system (x, y, r) . 
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IV. EXPERIMENTAL TECHNIQUES 


4.1 Limiting Streamline Angle Measurement 

The limiting streamline angles on a rotating body can be determined 
from the ammonia streak method described in Reference 36. A small 
amotmt of ammonia gas is released through the blade static pressure taps, 
while the rotor is in motion, and the trace of the gas is recorded on an 
ozalid paper. It is important to realize that, in order for the above 
procedure to be valid, the ammonia gas must remain in the collateral 
region, which means that the ammonia gas coming from the orifice on the 
blade surface must have a very low velocity. Typical traces of limiting 
streamlines derived by this method are given by Jabbari (37). 

Measurements of the limiting streamline angle are carried out at 
the predetermined tangential locations which coincided with the locations 
of the blade static pressure taps. At each tangential location, two 
measurements are made. Variation between the two measurements is found 
to be very small. The experimental accuracy of angle measurements is 
about 2° and for small cross flows which for example occur in the 
channel for open throttle case, the percentage error can be large ('\^10%) . 

4.2 Blade Static Pressure 

The flow in a channel differs from those around a single blade in 
two aspects: (a) The pressure gradient that exists in a channel changes 
the boundary layer behavior, local shear stress, limiting streamline 
angle, etc. (b) The interaction between the boundary layers on the two 
adjoining surfaces gives rise to complicated profiles arising out of 
mixing. In view of this it is important to have a knowledge of the 
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static pressure measureioents on the blade. Blade static pressure taps 

are connected to FTD by vinyl tubing and the static head measured by 

the micro-manometer. A centrifugal force correction term Is required 

to determine actual blade static pressure. 

If h Is the height of water column measured by the manometer, the 
m 

actual static head (h) Is given by 

2 

h = h^ + -- - (r - t /) (96) 

m p 2g o 

m o 

where p is the density of manometer fluid, r is the radius of 
m o 

rotating shaft (1") used in PTD, and r is the radius of the static 
pressure tap under consideration. 

The blade static pressure measurements were carried out at ten 
radii and at five chordwise locations, x/c - 0,1, 0,3, 0,5, 0.7, and 
0.9 on both leading and trailing surfaces of the blade, 

4, 3 Transition 

The boundary layer near the leading edge is laminar. The extent 
of the laminar region and the location of transition is determined by 
the "China Clay Technique", The blade surface is coated with a smooth, 
thin white film of china clay dissolved in acetone solution. The film 
is sprayed with ether by a spray gun and the rotor is brought to design 
speed. The evaporation takes place faster in the turbulent region than 
in the laminar region and is indicated by the whiteness of the former 
region compared to darkness of the latter region. The rotor is stopped 
after about three minutes and the laminar-turbulent transition line is 
obtained. Details of this method are described by Jabbarl (37) . The 



results are shown In Figure 16 for the case, when the rotor Is operated 
at large pressure gradients (with the throttle at a flow coefficient of 
0.05) . 


4.4 Skin Friction Coefficient 

One of the most Important problems In the three-dimensional boundary 
layer flow Investigations Is the determination of the wall shear stress. 
In two-dimensional boundary layer flows, the wall shear stress has been 
determined by various direct and Indirect methods. The freqiiently used 
direct methods Include the Preston tube method, the heated element method 
(HEM) and the floating element method (FEM) , 

Because of the presence of the centrifugal force, vibrations In 
rotation, use of FEM and HEM methods are severely limited In Its 
application for wall shear stress measurements on a rotating blade. 

The wall shear stress Is therefore evaluated by means of the 
Preston tube method. Validity of Its use for three-dimensional rotating 
turbulent boundary layer flows Is discussed below. 

Following Preston (38) when a total pressure pitot tube (whose 
diameter Is smaller than the extent of validity of the law of the wall) 


Is used, the wall shear stress can be related to the dynamic pressure 
AP recorded by the pitot tube through the following relationship. 


F (^) 
4p 


where d Is the diameter of the Preston tube. 

This shows that the wall shear stress can be determined by the 
local dynamic pressure, once the function F Is known. Experimental 
verification of the validity of this approach to stationary three- 
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dlmexislonal turbulent boimdary layer flows was carried out by Fierce 
and Krommenhoek (39). Results obtained by the Preston tube method were 
con^ared with those obtained by the floating element method which gave 
the direct measurement of the shear force, and good agreement was 
reported. 

In the present case of rotating three-dimensional turbulent boundary 
layers, the law of wall (u/u^ =* l/K^(log y^ + C)) Is experimentally shown 
to exist (Chapter V). Skin friction coefficients obtained from the 
Caluser Plot of the data In Chapter V agree very well with those obtained 
from Preston tube measurements and hence justify Its use. 

For the evaluation of the local total pressure, two different 
sizes of the pitot tubes (one of .028 inch outside diameter with 0.6 
diameter ratio, and the other .065 inch in outside diameter with .07 
diameter ratio) are used. Both are approximately 1.5 Inches in length. 

Measurements of the total pressure near the wall are carried out 
at five tangential locations, x/c = 0.1, 0.25, 0.5, 0.75, and 0.9 on 
both leading and trailing surfaces of the blade. At each tangential 
location measurements are carried out at four radial positions, r = 11, 
13, 15, and 17 inches. The tip of the pitot tube is placed flush to the 
blade surface approximately 1/8 inch upstream of the blade static 
pressure hole aligned In the direction of the measured limiting stream- 
line angle. The other end of the tube Is connected with a plastic tube 
to the pressure transfer device through the hole on the hub and the 
total pressure Is recorded by the manometer. Using the blade static 
pressure, the local dynamic pressure Is obtained. Following the 
calibration procedure given by Patel (40), the local dynamic pressure 
Is related to the wall shear stress. 
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All measurements at open throttle are carried out using both tubes, 
and the results are compared. They are shown in Figure 17, It is seen 
that the two probes give almost identical values Indicating that the 
diameter of the tube does not affect the calculated value of shear 
stress. All the subsequent measurements are therefore performed using 
a 0,065 inch diameter Preston tube, 

4.5 Mean Velocity Components 

Both three hole disk and the triple sensor hotwire rotating probes 
described in Chapter III are used for the measurement of mean velocity 
components. In the case of open throttle, only the disk probe is used. 
For the throttled case with large pressure gradients, a hotwire probe 
is used at all the measuring stations and the disk probe is used only 
at x/c = 0.9 to compare the two methods of measurement. 

In Case I (Open Throttle) the mean flow velocity profiles inside 
the rotor blade passage are obtained by measuring total pressure and 
the radial flow angle (a) with the three hole disk probe. Measuring 
grid is 5 X 5 X 40 in tangential, radial, and blade to blade directions 
respectively. 

The chordwlse velocity component u (parallel to the blade surface) 
and the radial velocity component w (cross flow) are: 



and 
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Figure 17. Radial and Chordwise Distribution of Skin Friction Coefficient (C^ on 
Leading and Trailing Surfaces of Flat Plate Inducer Described in Section 3.1 (Open^Throttle) . 



sin a 


(98) 


w 




where U. Is the blade speed at the tip and p Is the static pressure 
bt 

interpolated linearly from the blade static pressure measurements. 

For the throttled case, all the three mean velocity components u, 
V, and w are obtained by using the three sensor boundary layer hotwire 
probe described in Chapter III. The measured output voltages of the 
three channel hotwire anemometer, used in the hotwire setup (Figure 15) 
are solved numerically for the mean velocity components u, v, and w 
losing the appropriate calibration curves and the constitutive equations 
of the hotwires. The equations and the detailed method of solution is 
given by Gorton (23) . The method essentially relates the actual flow 
velocity components to the total effective cooling velocity on a finite 
heated wire in an arbitrary orientation of axes. The effective cooling 
velocity is further related to output voltage of the anemometer through 
the use of King's (41) Law. The equations of mean velocity components 
may be written in matrix form as; 

u. u. 

^ ^“in^(3,N) ^®J,ij^(3,3) ^ ^(3,N) 

(99) 

and 


iln^(3,N) 


= {e" J 


oV (3,N) 


{B^(n)}(i^3) 


^ ^(3,N) 


( 100 ) 


where, subscripts il, n, and i denote the wire, the observation point and 
the velocity component in the reference coordinate system respectively, 
u^^is the principal component of flow velocity, which in the present 
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are 


coordinate system is in the streanwise direction x. a„ , and 6« . . 

iti Xrij 

functions of direction cosines of the £-th wire which includes the 

correction factor for deviation from the Cosine Law (k, ) . E is the 

1 o 

output D.C. voltage at zero flow velocity, and are the 

measured values of voltage, slope of the hotwire calibration curve and 
effective mean cooling velocity for the Z wire and at nth observation 
point, u^ is the mean velocity component (e.g. u, v, w) , The numbers 
in parenthesis outside of a matrix represent its dimensions, e.g. 

jj) ® matrix with three rows and N columns, N is the total 
number of observation points in a measuring set. Indices 1, j, A can 
take values 1, 2, and 3. ^ 12 * °^13 Gorton’s terms 

®5* ®6* ^7 ^ 1 and 322 = ^33 " ^9 ^23 “ ^10’ 

other coefficients are zero. The values of functions . and 3« . . for 

xi xij 

the hotwire probe used in all measurements, are given in Appendix A. 

The computer program of Gorton (23), valid for a three sensor 
hotwire probe in a generalized non-orthogonal coordinate system is used. 
The program is modified to incorporate the following corrections to the 
hotwire data. The correction for deviation from the Cosine Law has 
already been included by Gorton (23) . 

1, Effect of changes in ambient fluid temperature during the course 

of the experiment and in calibration: The aiobient temperature drifts as 

much as lO^F in a measuring set and significantly affects the hotwire 

measurements. The hotwire equations are modified to Include this 

1 

correction as described in Appendix B. 

2, Effect of finite distances between three wires in a non-uniform 

flow: This effect is due to spatial resolution of the probe and can be 

quite significant in non-uniform flows, such as boundary layers and 
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wakes. Since the measurements have been done at more than one location 
in a traverse, this correction is incorporated by interpolating measure- 
ments from adjacent points, 

3. Effect of the variation of the slope of the calibration curve 
due to the deviation of hotwire equations from King’s Law; The hotwire 
output voltage does not follow King’s Law exactly i.e, the slope of the 
calibration curve is not constant with effective cooling velocity but 
varies according to the magnitude of cooling velocity. Since the cali- 
bration is performed over the entire range of effective cooling velocity 
anticipated in the measurements (0-80 ft/sec) this correction is 
incorporated by using the local slope, corresponding to the effective 
cooling velocity of a hotwire at any observation point, 

4. Effect of probe slip-ring contact resistance: It is small but 

can be incorporated if the anemometer output D.C. voltage is known at 

zero velocity, with and without the probe slip-ring contact resistance. 

2 

The hotwire constants of King’s Law are modified by a factor , 

where E^^ E^^ are D.C. output voltages of the anemometer at zero 

velocity, with and without the probe slip-ring contact resistance. 

5. Effect of the deterioration of the wires due to oxidation and 
prolonged use; This error can be made very small by operating the wire 
at optimum overheat ratio. An overheat ratio of 0.8 suggested by 
Wyngaard and Lumley (42) for velocity measurements is used. The wires 
are calibrated before and after each set. The variations in constants 
of the hotwire equations is found to be negligible. 

6. Effect of misalignment of the probe in the reference direction 

(x) : This error is due to inaccuracy in aligning the probe in the 

reference direction. For boundary layer flows, this error can be 
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minimized by using the well known fact that the velocity normal to the 
surface Is very small and goes to zero at the surface. The reference 
axis Is shifted by a small angle such that this velocity goes to zero 
for measurements nearest to the surface. 

A detailed discussion of these and other errors In hotwire 
measurements Is given In Appendix B. In addition the computation time 
was Improved by a factor of 10, by Improving the convergency rate and 
optimizing the calculation procedure. 


4.6 Turbxilence Intensity and Shear Stress Components 

All the six components of Reynolds stress tensor (u^*u^') at any 
point In the rotor passage In the rotating coordinate system (x, y, r) 
are obtained from the six measured values of the fluctuating voltage 
correlations three mean D.C. voltages (1, j 1 to 3), 

using the computer program of Gorton (23) . The fluctuating component 
of effective cooling velocity for ilth wire at nth data point may be 
written as: 


^^iln^(3,N) 


= + r 


n 


Jllj u 


^ ^in ^^^Jlli^(3 3) ^ 
In (3,3) (3,N) 


u ’. u’ 
in in 

'^In 


u; u' 
In in 


u 


In 
( 101 ) 


(3,N) 


The coefficients and nre same as occurring In Equation (99) for 

mean effective cooling velocity Q. The corresponding coefficients 
{r^lj} , for Jl = 1 In Gorton *8 form are given by 


r ssa r =sft r =r 

22 11* 33 12* 23 32 


a^^ and all other P's are zero. 


The coefficients of the first term can be symbolically expressed as. 
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u 


“jli u 


In 


(summation of j Is Implied) ■ 


( 102 ) 


Neglecting the last term in Equation (101) which is of second order, the 
matrix equation for Reynolds stress components Is: 




(103) 


where quantity (q^ q ) is the correlation between fluctuating cooling 
^ Jon mn 

velocities of Jlth and mth wires at nth observation point and Is related 
to output a.c. and d.c. voltages as follows: 




16 Eo E 
ion mn 




iln 


mn 




(104) 


The quantities B^j^(n) and etc. are respectively mean D.C. 

voltage, slope and mean effective cooling velocity of the ilth wire at 
the nth data point and N is the total number of data points in a 
measuring set. 

The matrix Equations (99) , (100), (103), and (104) are solved 
successively using the known values of voltages and calibration data, by 
using Newton-Raphson convergency scheme as described by Gorton (23) 
until the convergency criteria given below is reached. 

< 0.005 U, . (105) 


(P) _ „ (p-i) 
i X 


bt 


where u, and u. are values of a mean velocity component u. at 
1 i 1 

pth and (p-1) th iteration respectively and is the peripheral 
velocity of the rotor blade at the tip. 
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4.7 Spectrum of Turbulent Stress Correlations 


The spectrum of a Reynolds stress is defined as 

+00 

u^*Uj* “ / Co((jli) doi (106) 

•00 

where oi is the filtering frequency and Co(aj) is the amplitude of 
Reynolds stress at that frequency. It is therefore a measure of the 
contribution of various frequencies to the correlation and can also be 
used to check the contribution due to any spurious signal such as those 
caused by variation in main line frequency, probe variation, etc. The 
spectrum of all components of stress tensor u^’uj’ are obtained by 
using the Hewlett Packard 302A Wave Analyzer and Balantyne Mean Square 
Meter in the hotwire set up (Figure 11) , The spectrum in wave number 
space (k as 2 ttoj/U) is given by 

+00 

Uj^'u ’ = / iKk) dK (107) 

~ux> 

where (|)(k) is related to the measured Co(w) by 

4)(k) = (U/2it) Co(w) (108) 

U is the total mean velocity of the flow in the reference coordinate 
system (x, y, r) , 

The experimental value of Cospectrum at any frequency is obtained 
by filtering both u^' and u^ * at the same frequency, multiplying and 
averaging to give (Margolls (43)), 


u 'u ’ I = 4>..(w) dw 

^ J CO d 


(109) 


where (co) is the ar^>litude of the spectrum given in general by 

4>ij(w) » Co^j(o)) “ i Q^j((o), where i - / -1 (110) 
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is therefore the real part of the spectrum and the Imaginary part 

is called the quadrature spectrum. Since the time delay In u^' and 

u . * is zero in the present case, Q. . = 0, therefore, 
j 


<|)ij(aj) = Co^j(uj) = Measured Amplitude of <j) at O). (Ill) 

4)i^(oj) is calculated from the corresponding output voltages (e^e^) 
filtered and integrated at various frequencies, using equations (109) 
and (99) to (104) as described in Section 4.6 for the measurement of 
Reynolds stresses. The same computer program can therefore be used to 
obtain the spectrum of all turbulence quantities (u^’u^*). 

The spectral measurements are carried out only at midchord and at 
a radius of 17.25 inches. The main objective of the spectral measure- 
ments is to detect if there is any contribution to the electrical output 
signals from the spurious signals such as line frequency, blade passing 
frequency and its harmonics, natural frequency of the probe, etc. It 
is also required for the determination of the cutoff frequency (so as 
to eliminate contribution from the electrical noise above that 
frequency) of the filtering circuit in the hotwire set up (Figure 11) . 

It is found from the spectral measurements reported in Chapter V that 
there is no significant contribution to the hotwire output signals from 
the above mentioned spurious electrical signals. 
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V. EXPERIMENTAL RESULTS AND COMPARJISON WITH PREDICTIONS 


5.1 Data Without Throttle 

The experimental data Is first obtained by operating the rotor 
at the maximum flow rate, i.e. without usings a throttle at the exit. 

The flow in the rotor channel, therefore experiences only small pressure 
gradients and the major effects of rotation can be studied. The 
limiting streamline angle, blade static pressure and skin friction 
coefficient are measured both on leading and trailing blade surfaces. 

The mean velocity profiles are measured inside the rotor channel in a 
meas"ri .g grid of 5 x 5 x 40 in chordwise (x) , radial (r) and blade to 
blade (y) directions respectively. For this case, only the mean 
quantities which define the three-dimensional mean flow field have been 
measured and are described below. 

5.1.1 The limiting Streamline Angle . 

Radial variation of the limiting streamline angle is shown in 
Figure 18 for various cl* rdwise locations. It is seen that, except near 
the trailing edge, the limiting streamline angle on the leading blade 
surface behaves quite differently from that on the trailing blade 
surface. It is interesting to note that at x/c =* 0.1, the limiting 
streamline angles on the leading surface is nearly twice the value on 
the trailing surface. As x increases, this trend reverses itself, and 
at x/c = 0.5, the limiting streamline angle on the leading surface 
becomes smaller than that on the trailing surface. The flow is laminar 
for x/c < 0.1, where the limiting streamline angles are larger than the 
turbulent region (x/c > 0.1). The limiting streamline angles are 
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Figure 18. Radial and Chordwise Distribution of the Limiting 
Streamline Angle (Open Throttle) . 




generally hls^er than those for a single blade reported by 
Lakshmlnarayana et al. (17). 

The chordwlse variation of the limiting streamline angle (Figure 

18) shows that near the tip, the limiting streamline angle Increases as 

X Increases on the trailing surface while It takes a nearly constant 

value on the leading surface. For the radial distance less than 16 

Inches, the behavior of the limiting streamline angle Is different from 

that very near the tip. As x Increases, e Increases to reach a 

U) 

maximum value and then decreases on the trailing surface. This behavior 
Is quite different from the behavior of the limiting streamline angle 
for a single rotating helical blade which showed that the limiting 
streamline angle Increased monotonlcally with Increase in x. 

It must be pointed out here that for a single rotating helical 
blade the chordwlse pressure gradient Is zero, and the experimental 
Investigation of Lakshmlnarayana et al. (17) showed this to be the case. 
For the present investigation, blade static pressure measurements 
(Figure 19) have shown that there exist pressure gradients, throughout 
the passage. It Is Interesting to note that the behavior of the 
limiting streamline angle Is similar to the behavior of the blade static 
pressure. This can be seen when Figure 18 Is compared with Figure 19. 
This may be Interpreted as the Indication that the limiting streamline 
angle Is closely related to the blade static pressure distribution. 

Thus, In order to evaluate the limiting streamline angle It Is necessary 
that the pressure gradients be known.' 
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5.1.2 Blade Static Pressures 


The blade static pressure measurements are carried out at 11 radii 

and at five chordwlse locations, x/c « 0.1, 0.25, 0.5, 0.75, and 0.9 

both on leading and trailing surfaces of the blade. The results are 

plotted In Figure 19 for x/c « 0.1, 0.5, and 0.9. 

The radial variation of the blade static pressure coefficient at 

various chordwlse locations shows that ijj takes nearly a constant value 

s 

except near the tip where It changes considerably, especially near the 
trailing edge. 

The chordwlse variation of the blade static pressure coefficient 

shows that the static pressure distributions on the leading surface 

behaves quite differently compared to that on the trailing surface. 

While i|> takes the minimum value at the midchord location on the leading 
s 

surface. It takes the maximum value on the trailing surface. 

The existence of a finite loading on the blades can be explained 
on the basis of the fact that the boundary layers on the leading and 
trailing surfaces grow differentially giving rise to asymmetrical flow 
around the blade chord. 

5.1.3 Skin Friction Coefficient 

The skin friction coefficient Is obtained from blade shear stress 
measurements carried out at five radii (r - 11, 13, 15, 17, and 18 
Inches) and on both surfaces. The five chordwlse locations (x/c » 0.1, 
0.3, 0.5, 0.7, and 0.9) both on leading and trailing blade surfaces are 
shown in Figure 17. 

The skin friction coefficients derived from the velocity profile 
measurements reported In the next section, are plotted In Figure 20 In 
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Caluser Chart form l.e. u/u versus y"*" at x/c = 0.3 and 0.5. These 

e 

values are In close agreement with the values derived from the direct 

measurements shown in Figure 17. 

The skin friction coefficients are plotted as a function of = 

Reynolds number based on local momentum thickness 0^^ in Figure 21(a). 

The skin friction coefficient is found to decrease with Reynolds number 

Rq as a power law of index (-0.268). Effect of the limiting streamline 
®11 

angle (e ) and the shape factor (H) is best represented by the 
w 

expression 

c. a (10)”^*^^®^ (1+0.52 /e (x-x ) /c ) 

f w c 

as shown in Figure 21. Thus the skin friction coefficient correlation 
for rotating turbulent boundary layers with mild pressure gradients are 
best represented by the following formula 

c^ = 0.172 (R. )"0*268 ^q^-0.678H ( 2 +q 52 /e (x-x ) /c ) 

f 6 ^^ w c 

( 112 ) 


where 

c^ = T /(l/2)pu ^ . 
r o e 

It should be noted that this correlation is similar to that of a single 

rotating disk of Lakshmlnarayana et al (17) , where the effect of 

1/2 

rotation parameter (£ x/c) is discussed. The term in the bracket 

w 

in Equation (112) represents the effect of three-dimensionality in the 
flow. It should be noted that this correlation is similar to Ludwig and 
Tlllmann^s (28) correlation for a stationary turbulent boundary layer 
with pressure gradients l.e. 


c. = 0.246 
f 


( 10 ) 


-0.678H 


(Rf) ) 
®ii 


-0.268 


(113) 
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Equation (112) has a different constant and an additional parameter 

(1 + 0.52 /e (x-x ) /c ). However, It should be noted that since Is 
w t 11 

found to decrease appreciably In the rotating case due to an Increased 

radial migration of the flow, Is much less for the same free stream 

11 

velocity and the pressure gradients and therefore the effect of 
rotation la to Increase the skin friction coefficient through a decrease 
of compared to the equivalent stationary case. 

®u 

5.1.4 Mainstream Velocity Profile . 

Figures 22-26 show the distribution of chordwlse velocity 
component 'u' across the passage at various radial and chordwlse 
locations. From the profile distribution at x/c ** 0.1, it is evident 
that the flow has already become turbulent over most of the radial 
locations. The boundary layer is thicker on the trailing side (as much 
as one Inch near the tip). This may be due to inlet flow incidence. 

The boundary layer is thin and laminar near the hub on the leading 
surface (at the hub, the Reynolds number based on the peripheral 
velocity, J^r and distance from leading edge is approximately 0.8 x 10^), 
The profile distribution at subsequent downstream measuring 
stations, (Figures 23 - 26) shows that the boundary layer Is turbulent 
at all the radial locations on both the surfaces and grows rapidly 
downstream. It should be remarked that the free stream velocity, unless 
the two boundary layers have merged is found to be the same as the blade 
speed (^r) . 

The boundary layers at the tip grow rapidly due to radial flow 
migration. Fully developed flow occurs at x/c “ 0.30 for R ■ 0.98 and 
this region extends Inward as the flow travels downstream. One of the 
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Figure 22. Blade to Blade Distribution of Streawlse 
Velocity Component at x/c = 0.1 (Open Throttle). 
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Figure 24. Blade to Blade Distribution of Streanwise Velocity 

Component at x/c 0.5. 
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Figure 25. Blade to Blade Distribution of Streamwlse 
Velocity Conponent at x/c -0.7 (Open Throttle). 
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Figure 26. Blade to Blade Distribution of Streamvlse 
Velocity Component at x/c » 0.9 (Open Throttle). 
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unique features of the channel flow Is the Interference region, where 
the boundary layer profiles are unconventional. This can be discerned 
for x/c > 0.3, near the tip (Figure 26). A wake type of profile exists 
near the mldpassage. As explained earlier In Chapter II, this Is brought 
about mainly by the Interaction of blade boundary layers and the annulus 
wall. The resulting Inward flow near the mldpassage (Figure 2) produces 
this type of profile. The Invlscld flow region disappears and therefore 
no definite boundary layer characteristics can be obtained In terms of 
boundary layer thickness parameters (6, 6^*, etc.). It Is 

obvious that in this region the momentum Integral technique can only 
give a qualitative picture of the flow. 

The boundary layer Is generally thicker near the trailing surface 
except near the trailing edge (Figure 26). 

5.1.5 Comparison with Power Law Profile . 

In Figure 27, the mainstream velocity is plotted against the 

nondlmenslonal distance from the blade wall (y/6). Conq)arlng this with 

1/n 

the Power Law profile of appropriate Index (u/U^ “ (y/5) ) it appears 

that the experimental data fits this representation reasonably well. 

The agreement Is poor near the tip because the boundary layer here Is 
affected by the tip clearance and the annulus wall effects. The value 
of exponent varies from 1/6 to 1/11, except near the hub at x/c “ 0.3 
where the values are 1/4 and 1/5. Here the botmdary layer appears to be 
In transition (Reynolds number based on peripheral velocity and 
distance from the leading edge Is 3 x 10^) . 

Comparison with Cole*s Profile: In Figure 28, the mainstream 

velocity component u Is plotted In the form of Cole's Law, 
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Figure 27, Power Law Representation 


Velocity Component (Open Throttle) 


die 







(114) 


u/u^ = l/K^(log(yu^/v)+Tiuj(y/6)) + C 

where is the friction velocity, is von Karman constant, U) is a 
wake function, 1T is a parameter to account for the effect of Imposed 
pressure gradients, and C is a constant. 

It is seen from the plots that there does exist a law of the wall 

+ 

region (u/u^ = 1/k log yu^/v + C) up to y < 200-300 depending upon the 
nature of the blade surface (trailing or leading) and the location of 
the measuring station. The mean value of constants in the law of the 
wall vary from 0.40 - 0.42 for the slope constant K and 4.5 to 5.5 for 
the constant C. The wake component 03(y/6) of the Cole's Law has a 
different form than the sinusoidal form for stationary two-dimensional 
flows. This may be due to the fact that the entrainment process, which 
is dominant in the outer part of the boundary layer, is different for 
the boundary layer on an enclosed rotating surface than for a stationary 
two-dimensional surface. Value of ir which is characteristic of the 
Imposed pressure gradients is found to vary from -0.25 to 0.55. It is 
positive on the trailing side and has a small negative value on the 
leading side of the blade. It decreases with x up to x/c < 0.5 and 
later Increases slowly. This is in accordance with the trend of the 
observed blade static pressure distribution, as shown in Figure 19. The 
value of 7T increases with radius continuously in accordance with the 
radial pressure gradient. 

Distribution of Passage Averaged Chordwise Velocity; In Figure 29 
is shown the radial variation of mainstream velocity component averaged 
across the channel. The passage averaged value of u is found to be 
lower than the free stream velocity that exists in an inviscld flow 
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Streamwlse Velocity CoBqjonent (Open Throttle) . 
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Figure 30. Distribution of Shape Factor In the Rotating Channel 

(Open Throttle) . 





(u =» fir) . The difference (fir - u) Increases as the flow proceeds 
further downstream of the channel, especially after x/c = 0.5. The 
largest difference Is near the tip where the two surface boxmdary layers 
have interacted with the annulus wall boundary layer, and for x/c >0.5 
it extends up to mid radius (R 'v 0.8). 

This large reduction in mainstream velocity near the tip is due to 
the mutual interference of the boundary layers on the two surfaces with 
the annulus wall boundary layer, which are now completely submerged into 
each other and is accompanied by radial inward flow in this region. 

This is in conformity with the increase in radial pressure gradients 
observed near the tip region. The value of shear stress is not negligible 
near the center of the passage and the turbulence level in this region 
is maintained at the expense of mainstream velocity gradient in the 
normal (y) direction. These results are confirmed by turbulent stress 
measurements reported in Section 5.3.5 

Shape Factor (H) : The chordwise distribution of shape factor H 

obtained from the measurements is plotted in Figure 30. The value of 
H is higher on the leading side (1.30 to 1.50) conq>ared to that on the 
trailing side of the channel (1.10 to 1.30). It is nearly constant near 
the mid radius but changes rapidly at the hub and the tip due to the 
interference effects of the hub or annulus wall boundary layers. The 
value is more at higher radii due to larger radial pressure gradient. 

It should be remarked here that the measured values of H are only 
qualitative and are liable to some error (due to interpolation of 
velocity profiles near the blade surface to obtain H) particularly 
where the boundary layer is thin. Since this set of experimental data 
has been obtained by running the four-b laded channel at open throttle 
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and hence with mild pressure gradients, the above results show that It 
Is valid to take H constant In x direction In predicting 0^^ and In 
this case. 

5.1.6 Radial Velocity Profile 

Figures 31 and 32 show the observed radial velocity dlstrlbut(lon 
across the blade passage. The radial velocity Increases as the wall Is 
approached. It should decrease very close to the wall, but since the 
measurements could not be taken for values of y^ < 100 (due to finite 
size of the probe) , only the outer region of the boundary layer flow has 
been scanned. Near the centerline of the channel, the radial velocity 
Is small (S^ 1 ft/sec) and Is within the experimental error (± 1 ft/sec) , 
except near the tip. The radial Inward flow Is observed near the tip 
for x/c > 0.3. The radial outward flow region Is confined to blade 
boundary layers, being larger near the trailing surface for x/c < 0.5. 
For x/c > 0.5, the leading surface boundary layer has larger radial 
outward component. It Is clear that the radial flows are appreciable 
In the rotating blade channel and are higher than those of a single 
rotating blade reported by Lakshmlnarayana et al. (17). 

The observed radial velocity component w Is plotted In Figure 33 
as streamline parameter e (local radial to mainstream velocity ratio) 
normalized by the local measured streamline angle > l*c. z/z^y 
versus nondlmenslonal distance from the wall (y/6) , Measurements near 
the outer edge of the boundary layer where the radial velocities are of 
the same order as the experimental error 1 ft/sec) are not shown for 
clarity. A comparison with Mager's profile (Equation 74) Indicates that 
the experimental data fits this form quite well except near the hub or 
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tip. In the interference region the value of e/e is higher than Mager's 

w 

profile and near the hub the values are, in general, lower. The 
discrepancy is due to the effects of endwalls which interfere with the 
blade boundary layer producing complex flows. In the tip region, where 
the radial Inward flows are present, Mager's representation fails 
because of the reversal of radial velocity profile (s type of profile) 
Inside the boundary layer. This occurs at x/c » 0.5, and R > 0.93, 
x/c = 0.7, and R > 0.86 and at x/c =0.9 and R > 0.86. 

Plotted in Figures 34 are the hodograph plots of velocities on the 
leading and trailing surfaces respectively. The radial velocity measure- 
ments are inaccurate near the hub for x/c = 0.3, since the boundary 
layer thickness is very small there and therefore this data has been 
omitted in this polar plot. The data is compared to following predicted 
forms for the outer layer. 

(i) Johnston's Triangular Model w/U^ = ^ 2 ^^ ~ (115a) 

(ii) Author's Model w/U = C- /(2/k-sin3) sin3/RG-(u/U ) (115b) 

e i ® 

It is clear from the figure that Johnston’s straightline fit 115a doesn't 
agree with the experimental data (except at some locations near the hub, 
which is coincidental) . The agreement with Equation (115b) is quite good 
for most of the data, if the value of is allowed to vary with x. In 
the region where the two blade boundary layers have not merged i.e. for 
x/c < 0.5, it has been experimentally shown in Figures (23) and (24) 
that the velocity at the edge of boundary layer (U^) ^^r and therefore 

k = RG - 1.0 and 

w/U = 0.2 /(2-sin3)sln3 - (u/U )^ 
e ® 

gives the best fit. 
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Equation (31) as well as the experimental data Indicate that It Is 

not necessary for the radial velocity w In the rotating boundary layers 

In a channel to go to zero when u » U and BU /dy 0 as has been 

e e 

observed In the case of boundary layers In the stationary passages. 

Polar plot of the velocity profiles near the blade tip are shofwn In 
Figure 35. Since the radial velocities are comparatively small at x/c 
■ 0.7, no definite conclusions can be drawn with regard to the profile 
shape. However, at x/c « 0.9, (Interference region - Figure 2), the 
radial Inward velocities as well as the defect In mainstream velocity Is 
appreciable. The data Is compared with the Equations (115a) and (115b) 
with a negative constant In the "Interference region". It Is evident 
from the figure that Equation (115b) gives the best fit except near the 
center of the channel section. The straight line form (115a) Is adequate 
If C 2 Is assumed different on the two surfaces. 

The existence of such a phenomenon In turbomachinery has not been 
reported before. It should be remarked here that the mixing losses 
associated with the Interference region has never been taken In to 1 
account In developing end loss correlations in turbomachinery. This 
accounts for some of the discrepancy observed between endwall 
correlations (based on cascade data) and the endwall losses measured 
at the exLt of turbomachinery rotors. 

5.1.7 Gomparlson of Experimental Results with Theory 

The momentum Integral technique described In Chapter II has been 
used to predict the momentum thickness (^^^-j^) sud limiting streamline 
at various r and x locations for the four bladed flat plate Inducer. The 
measured static pressure distribution Is used to obtain RG, T6 and k. 
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The predicted variation of and 0^^ for the present case as well 
as single blade case are also shown In Figures 18 and 36 respectively. 
For the present case of the flat plate Inducer run at open throttle » the 
values of 6^^ obtained from the experiments at first Increase rapidly 
up to x/c « 0.3 and decrease or Increase slowly further downstream of 
the channel depending upon the Imposed pressure gradients. The values 
are higher on the trailing surface compared to the leading surface which 
Is due to higher pressure gradients (Figure 19) found on the latter 
surface. 6^^ Increases with radius for most of the channel » except near 
the trailing edge of the trailing blade surface. The predictions agree 
well with the measured values except near the leading edge, where the 
extent of laminar boundary layer growth Is assumed to be known. The 
location of transition Is determined experimentally. Compared to the 
boundary layer growth on a single blade rotating In stationary fluid, 
the boundary layer growth Is much smaller In a channel due to the 
presence of radial and chordwlse pressure gradients. The difference 
between the observed and predicted boundary layer growth Is partly due 
to an error In estimating the extent of Initial laminar boundary layer 
and partly due to an error In estimating the momentum thickness 9^^ from 
measurements. A study of Equation (70) for 0^^ shows that the effect of 
adverse chordwlse pressure gradients Is to Increase the boundary layer 
thickness, while the effect of positive radial pressure gradients Is to 
decrease It compared to the single blade case plotted In Figure 36. It 
can therefore be concluded that In the present case, the effect of 
radial pressure gradient Is much larger than the chordwlse pressure 
gradient. This Is confirmed by a large decrease In 0^^ near tip radius 
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where the radial pressure gradients are found to be much larger than 
the chordwlse pressure gradients. 

The measured and predicted values of the limiting streamline 
angle are compared in Figure 18. In the turbulent region i.e, x/c ^ 0.1 
the values Increase or decrease with x in accordance with the pressure 
gradients. The value of decreases with Increase in radius and 
reaches a mlniiium at the tip. The magnitude of is higher on the 
leading side indicating that the boundary layer is more skewed and that 
the cross flow is larger compared to that on the trailing side of the 
channel. The predictions agree reasonably well except near transition. 
The value of in the laminar region is obtained from the approximate 
solution of Banks and Gadd (7). The numerical results are subject to 
some error because the transition is assumed to be sudden and there is 
an order of magnitude of difference between the values of limiting 
streamline angle laminar and turbulent boundary layers 

respectively. An important difference between the single blade results 
of Lakshminarayana et al. (17) and the rotating channel is that the 
values are higher in a rotating channel and there is a radial variation 
due to the presence of varying radial pressure gradient. 

An examination of Equation (71) indicates that is very sensitive 
to pressure gradients and its value Increases in the presence of positive 
radial and chordwlse pressure gradients. This accounts for larger values 
of measured in a channel as compared to a single blade, plotted in 
Figure 18. The extent of radial flows and boundary layer growth is very 
sensitive to radial and chordwlse pressure gradients. This conclusion 
Is confirmed by the comparison of boundary layer growth on a single 
blade and in a channel shown in Figure 36. 
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5. 2 Effect of Throttle and Exit Flow Measurements 

In order to Investigate the overall performance of the four b laded 
flat plate Inducer at varlotis flow coefficients, a throttle is installed 
at the exit as shown in the test facility (Figure 8) , The second 
objective is to determine the operating conditions for conducting the 
measurements Inside the rotor so that the pressure gradients in its 
blade passages are of the same order of magnitude as those found in the 
actual rocket p\amp inducers (e.g. Lakshmlnarayana (21, 45)). 

The throttle clearance is varied and annulus wall static pressure 
measurements are carried out from inlet to far downstream. The flow 
rate, stagnation pressure and flow angle are carried out 14 inches 
downstream of the rotor. 

Figure 37 shows the static pressure distribution measured on the 
annulus wall with different values of throttle clearance. It is seen 
that the annulus wall static pressure increases as the throttle 
clearance decreases. The negative pressure at the leading edge of the 
blade is considered to be caused by the blade blockage effect. Increase 
of the blade static pressure at the trailing edge is considered to be 
caused by the increase of flow area due to the absence of the blade 
thickness. It is seen that the static pressure is constant downstream 
of the rotor assembly. It is Interesting to note that the axial pressure 
gradient on the annulus wall decreases as the throttle clearance 
decreases. The stagnation pressure rise coefficients at various throttle 
positions are plotted in Figure 38. Even at large throttle openings, 
appreciable pressure rise can be seen especially at the tip. This is 
caused by considerable boundary layer growth on | the blades, and the 
shear pumping effect. The outlet angles, measured with respect to 
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Figure 38, Distribution of Stagnation Pressure Coefficient 

(i/^^) at the Inducer Exit. 
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tangential directions are derived using a cobra probe. The flow 
reversal near the hub, occurs at flow coefficient of 0.0592 or less as 
indicated by angle measurements shown in Figure 39. The performance 
of the flat plate inducer, plotted in the form of - cp curve is shown 
in Figure 40, where is mass averaged stagnation pressure coefficient 
and (j) is passage averaged flow coefficient. 

For experiments at large pressure gradients, a throttle clearance 
of 0.7 inch corresponding to a flow coefficient of 0.05 and an average 

i 

stagnation pressure coefficient - 0.25, is chosen.' Measurement of 

stagnation pressure and flow angle are performed at stations 3 and 6 

7/16" and 14" downstream of the rotor trailing edge. Static pressure 

is measured both at the hub and annulus walls. In Figure 41, ip and 

X s 

distribution is compared to that of a four bladed cambered rocket pump 
inducer tested by Lakshminarayana (21, 45). It is seen that the two 
rotors have similar characteristics even though the value of for the 
flat plate rotor is about 25% lower. Axial and radial velocity profile 
distribution (Figure 42) show similar trends. However, the backflow 
region for the flat plate inducer is more because it is being operated 
at lower flow coefficients, 

i 

5. 3 Experimental Data with Throttle 

Actual turbomachinery rotors have large pressure gradients at 
design flow conditions. In order to simulate the same flow conditions 
in the experimental rotor, the flow is partially throttled at the exit. 

A throttle clearance of 0.70 inch corresponding to a flow coefficient 
of 0.05 is chosen so that its stagnation pressure distribution (operating 
characteristics) corresponds to an actual rocket pump inducer operating 
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Figure 41. Comparison of Radial Distribution of Stagnation Pressure 
Coefficient and Static Pressure Coefficient for Flat 

Plate and Cambered Bladed Inducers. 
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at its design conditions. All the mean flow quantities obtained for 
rotor operating without throttle are also measured in this case. In 
addition turbulent intensities and stresses and their spectra are 
measured to study the effects of rotation on the local turbulent 
structure of the flow in the rotor channels. A three sensor hotwire 
probe described earlier In Chapter III is used for all mean velocity 
and turbulence measurements in a measuring grid of 5 x 5 x 39 in radial, 
chordwise and blade to blade directions respectively, (except at x/c = 
0.1, where measurements are carried out at 3 radii). 

Mean velocity profiles are also measured at x/c = 0.9 using the 
three hole disk probe (Figure 12) to compare the obtained data by two 
separate methods of measurement. The agreement between two measurements 
is found to be good. The experimental results of all measurements 
inside the partially throttled flat plate inducer are described below. 

5,3.1 The Limiting Streamline Angle 

Figure 43 shows both the radial and the chordwise distribution of 

the limiting streamline angle both on leading and trailing blade surfaces. 

The radial distribution shows a continuous decrease of e with radius 

w 

at all chordwise locations, except near the hub where the interference 
of the hub wall boundary layer with blade boundary layer is appreciable. 

At x/c = 0.1, the flow is laminar near the hub on the leading side and 
consequently results in large values of the limiting streamline angle. 

The chordwise variation of the limiting streamline angle (Figure 43) 
shows that in general the magnitude increases with x due to adverse 

chordwise pressure gradients. However, for values of x/c > 0.75, the 
limiting streamline angle increases on the leading surface and decreases 
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on the trailing surface In accordance with the pressure gradients : 
(Figure 44). The magnitude of limiting streamline angles are In general 
much higher compared to both single blade and open throttle data. This 
Is becaiise the pressure gradients (Figure 44) are much larger In the 
partially throttled case compared to the open throttle case (Figure 19) . 
cThe values are comparable to those for actual four b laded cambered 
rocket pump Inducers where the pressure gradients are also of the same 
order. This confirms an earlier conclusion, based on open throttle data, 
that the magnitude of the limiting streamline angle is directly affected 
by the lnq)osed external pressure gradients on the turbulent boundary 
layers In a rotating passage. 

5.3.2 Blade Static Pressures 

The radial distribution of blade static pressure at various 
chordwlse locations Is shown In Figure 44. The measuring points are 
the same as In the open throttle case. Unlike the open throttle case 
(Figure 19) , the radial variation of Is large. Its gradient also 
Increases continuously with the radius to a maximum value at the tip. 
Except In the entrance region (x/c = 0.1), the static pressure on the 
leading surface Is higher compared to that on the trailing surface. At 
the entrance, since the Inducer Is being operated at large Incidence 
(4° at hub to 2 ° at tip) , the leading side has a larger negative value 
of static pressure than on the trailing side of the blade surface. 

It Is also clear from Figure 44 that the static pressure rises 
rapidly with x up to quarter chord. Further downstream, for x/c > 0.25, 
the static pressure coefficient remains relatively flat on the 

leading surface even up to ninety percent of the chord. On the trailing 
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side ip increases continuously up to x/c - 0.75 and decreases there 

onwards. The pressure loading on the surface occurs mainly between 

x/c *■ 0.1 and 0.5, and is more pronounced at hub and tip radii. Very 

near the tip (R > 0.97), ip decreases at x/c = 0.1 on leading side due 

s 

to the presence of leakage and secondary flows at the tip. Due to 
large chordwise pressure gradients and incidence, the relative flow 
near the hub is found to be separated near the leading edge (e.g. at 
x/c = 0.1). 

5.3.3 Mean Velocity Profile 

Figures 45 to 49 show the distribution of chordv/lse and radial flow 
velocity components at various radial and chordwise ]acat:ions. Their 
values are obtained by using three sensor rotating hot,v7ira probe and 
hotwire equations as described in Chapter IV. The velocity cotoponent 
normal' to the blade is found to be very small throughout the channel 
(< 1 ft/sec) and has not been shown in these figures for the sake of 
clarity. In order to check the validity and accuracy of hotwire 
measurements in rotating system, these mean velocity measurements were 
compared to disk probe measurements at x/c = 0.9 at fi\'e radii. The 
measurements by both methods gave almost Identical valxiea of mean 
velocity profiles within the range of experimental error and hence 
establish the validity of the hotwire measurements in rotating frame 
of reference. 

Chordwise Velocity Profile: An examination of the stream^ise 

velocity at x/c - 0.1 shows that the flow has already become turbulent 
over most of the radii. The boundary layer is thin on the leading side 
and thick on the trailing side. The velocity profiles clearly show that 
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the flow is asymmetrical and that large velocity gradients exist across the 
channel. This is caused by large chordwise and radial pressure gradients 
that exist near the leading edge (as shown in Figure 44) in the present 
case. The inducer flow has been throttled to give a large pressure rise 
at low flow coefficient as shown in the performance curve of Figure 40, 

The profile distribution at subsequent downstream measuring stations, 
Figures 45 to 49, show that the boundary layer is turbulent at all radii 
on both the surfaces and grows rapidly downstream. It should be noted 
that unlike the open throttle case discussed in Section 5.1, there is a 
large diffusion of relative flow due to large positive chordwise pressure 
gradients (Figure 44) in the channel in the entrance region up to x/c =* 

0.5. Therefore, the streamwise velocity in the Inviscid region external 
to the blade boundary layers is smaller than the blade speed (^2r) . 

The boundary layers at tip grow rapidly due to radial flow migration 
and large unfavorable chordwise pressure gradients. Fully developed 
flow occurs at x/c *» 0.30 for R ^ 0.94. This region extends inward 
rapidly. In the interference region (defined in Figure 2), the boundary 
layer profiles are unconventional, as observed and explained earlier in 
the open throttle case in Section 5.1. As a result, the boundary layer 
thickness cannot be defined in a conventional manner and the momentum 
Integral analysis which assumed the monotonically increasing powerlaw 
profile for the streamwise velocity is inadequate. It gives only a 
qualitative picture, if the boundary layer thickness is defined as the 
distance from the surface, where the streamwise velocity component 
reaches a maximtua. 

An examination of velocity profiles shows that the boundary layer 
on the trailing surface is first thicker compared to leading surface up 
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to midchord. Further downstream, the trend is reversed and at x/c ■ 

0.5, the boundary layer is thicker on the leading side. This behavior 
is in accordance with the chordwise pressure gradients (Figure 44) . The 
pressure gradients on the leading surface is almost constant after 
midchord but it has a negative value which results in the acceleration 
of the mean flow on the trailing surface. The boundary layers on blades 
are in general thicker in the throttled case compared to those in the 
open throttle case. 

The mainstream velocity is plotted in Figure 50 against 
non-dimensional distance from the blade wall (y/6) at x/c *> 0,3, 0.5, 
and 0.7 and at radii R = 0.71, 0.86 and 0.93, both for leading and 
trailing side boundary layers. Measurements at hub and tip radii have 
not been plotted in this form, because the boundary layer behavior is 
unconventional due to end wall effects. As noted in the open throttle 
case, the measurements fit the power law profile representation well. 

The value of power index varies from 1/6 to 1/11 depending upon the 
nature of pressure gradients (Figure 44) , the value of power index 
decreases downstream from x/c = 0.3 to 0.7 particularly on the leading 
side, which has favorable pressure gradients in this region. 

Coles form of streanwise velocity has not been plotted because 
the measurements in the open throttle case plotted in this form have 
clearly shown that the usual sinusoidal type of wake profile observed in 
the outer region of stationary turbulent boundary layers is not observed 
in rotating turbulent boundary layer data. 

In Figure 51, the radial distribution of streankrlse velocity at the 
edge of the bomdary layers is plotted, both for leading and trailing 
side blade boundary layers, at five chordwise locations. Unlike the 
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Figure 50. Power Law Kepresentatlon of Velocity Profiles in the Blade Boundary Layer (Throttled Case) 
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open throttle case, there is a large decrease in the boundary layer 
edge velocity even in the entrance region. This is due to the fact 
that the Inducer in the present case is operated at a lower flow 
coefficient (0.05) compared to the open throttle case (0.065), which 
results in an appreciable pressure rise as shown in the performance 
curve of Figure 40. The difference is larger near the hub and the tip, 
where the hub and the annulus walls interfere with the flow. Since 
the inducer is operating at a large incidence (4^^) , the relative flow 
is separated at the leading edge and causing the observed reduction 
in the strearawise velocity. At the tip, the interference with the 
annulus wall produces radial inward flows and a reduction in the 
strearawise velocity as described earlier for the open throttle case. 

This is accompanied by an increase in radial pressure gradients at 
the tip as shown in Figure 44. The decrease is larger on the leading 
surface compared to the trailing surface, particularly in the entrance 
region, where the leading side encounters larger unfavorable pressure 
gradients. 

Radial Velocity Profile: The blade to blade distribution of the 

radial velocity profiles at different radial and chordwise locations 
are also shown in Figures 45 to 49. The radial velocities show the same 
behavior as noted in the open throttle case i.e. an increase in 
magnitude on approach to the blade surface. The magnitude of radial 
flows increase downstream reaching a maximum value near the trailing 
edge (x/c = 0.9). The radial velocities are first larger in magnitude 
on the trailing surface compared to leading surface but after midchord 
the trend reverses itself. Similar behavior is observed in the open 
throttle case and is due to the fact that the chordwise pressure 
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gradient on the trailing side Is favorable to the flow after the 
mldchord, while on the leading surface, the static pressure Is almost a 
constant. The magnitude of radial flows are larger as compared to both 
the single blade and the open throttle cases. This can be also 
discerned by conqparlson of the plot of limiting streamline angles in 
Figures 18 and 43. The values in the throttled case are much higher 
than the other two cases. This is mainly brought about by the presence 
of larger radial and chordwlse pressure gradients In the present case 
compared to zero pressure gradients In the case of the single blade and 

mild chordwlse pressure gradients in the open throttle case. 

% 

The radial velocities are plotted In a hodograph plane in Figure 

52, both on the leading and the trailing surfaces and are compared to 

the theoretically predicted form (Equation 115b) . The agreement is 

quite good at most of the locations. It is seen clearly from the data, 

that the radial velocity does not go to zero at the location where 

u = U . As explained earlier in section 21. , this is due to the fact 
e 

that centrifugal force is larger than the radial pressure gradient 
~ 2 

force Vq /r) in the rotor channel and results in outward radial flows 

even if 9u /3y = 0, at places where U sin$^f^r. 
e e 

5.3.4 Turbulent Intensities 

Blade to blade distribution of all the three turbulent Intensity 
components at various radial and chordwlse locations Inside the rotor 
channel are shown in Figures 53 to 73. The turbulent intensity 
components have been normalized by local total relative flow velocity (U) . 
Figures 53 to 56 show the measurements near the hub. Figures 57 to 65 
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Blade to Blade Distribution of Turbulent Intensity 
at R = 0,52 and x/c = 0.3. 
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Figure 56. Blade to Blade Distribution of Turbulent Intensity and Stress Correlations 

at R = 0.564 and x/c = 0.9. 
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Figure 60. Blade to Blade Distribution of Turbulent Intensity and Stress Correlations 

at R = 0.712 and x/c = 0.7. 
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Figure 65. Blade to Blade Distribution of Turbulent Intei 

at R = 0.866 and x/c * 0.9. 
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Figure 66. Blade to Blade Distribution of Turbulent Intensity and Stress Correlations 

at R = 0.911 and x/c = 0,1. 
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Figure 69* Blade to Blade Distribution of Turbulent Intensity and Stress Correlations 

at R «= 0.938 and x/c « 0.7. 
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Figure 72. Blade to Blade Distribution of Turbulent Intensity and Stress Correlations 

at R s 0.983 and x/c =0.7. 
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Figure 73. Blade to Blade Distribution of Turbulent Intensity and Stress Correlations 

at R “ 0.98 and x/c » 0.9. 



0 

Leading 

Surface 


show the measurements near the -midradius and Figures 66 to 73 present 
the data near the tip radius. Experimental accuracy In measurements Is 
estimated to be within 5% for all Intensity components. The analysis 
of errors In turbulence Intensities is described In Appendix B (section 
B.4). 

The magnitude of turbulence Intensities Is higher near the hub 
and tip regions due to the effect of end walls. The values are 
highest near the hub ('^^16% at x/c = 0.3 as shown in Figure 53) which 
operates at a maximum flow Incidence (4*^). The mean velocity profiles 
(Figs. 45 and 51) have shown the presence of hub wall boundary layer and 
this may account for the large values of turbulence intensities in this 
region. 

The free stream turbulence level at all radii in the entrance 
region is quite high. It varies from 20% near the hub to 6% near the 
tip at x/c = 0.1. The free stream turbulence level is found to decrease 
with distance inside the rotor channel. It decreases from about 8% 
(Figure 55) and 4.5% near the tip at x/c =0.7 (Figure 72). On the 
other hand, near the mid radius, the free stream turbulence level is 
almost constant at R = 0.7 (Figures 57 and 61). It increases from 3.5% 
at x/c =0.3 (Figure 62) to 5% (Figure 65) at x/c = 0.9 for R - 0.86. 
This is brought about mainly by the turbulent mixing of the flow in the 
channel in this region. 

The turbulent intensities increase on approach to the blade 
surfaces. This Increase is due to the production of turbulence by 
mean velocity gradients in the blade boundary layers . They should 
decrease very near the surface, in the viscous layer but since the 
measurements could not be taken for distances less than y*^ < 50 (due to 
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finite size of the probe) this phenomenon has not been observed 
experimentally. 

Except at x/c » 0.1, the magnitude of radial intensity component 

. 2 '2 2 

(w* ) is larger than both the streamwlsc (u' ) and blade normal (v' ) 
intensity components. This observed. behavior is opposite to known 

"^2 2 

behavior for stationary two-dimensional boundary layers where u* > w* > 
2 

v' and clearly indicates that the turbulence characteristics are 

markedly changed by the effect by the effect of rotation. Experimental 

data of Gorton (23) inside the rotor passage of a three bladed rocket 

pump Inducer and of Raj (46) In rotor wakes confirm this conclusion. 

Qualitatively, it is shown In Appendix C that, in compressors and 

pumps where the direction of rotation is opposite to that of the 

relative flow, the rotation yields an extra production term in the 

transport equation (Appendix C, Equation C-3) of radial intensity 
2 

component (w* ) and a suppression term in streatowlse direction provided 
u’w' correlation is negative. Shear stress measurements in Figures 53 
to 73 show that this is generally true inside^ the rotating blade 
boundary layers in axial flow machinery. 

2 2 

However, at x/c = 0.1, Figures 57 and 66 show that u* > w* > v' 
which is similar to the observed behavior in statlonairy two-dimensional 

boundary layers. An examination of transport equations show that it is 

2 ~'~2 

possible to have suppression of w* and production of u* if the 
correlation (u*w*) is positive. The turbulent stress data in Figures 
57 and 66 indicates that this is true at x/c - 0.1. 

Turbulence Intensity measurements shown in Figures 58 to 65 
represent data free of endwall and leading edge effects and is 
representative of flows in axial-flow rotors operating at design flow 

\ ■■ 
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conditions. It is seen that turbulence intensity components increase 
downstream with distance from x/c =0.3 to x/c = 0.9 on the leading 
side and decrease on the trailing side. This is mainly brought about 
by the imposed chordwise pressure gradients (Figure 44) in the passage 
in this region. These pressure gradients are very small on the leading 
side and large and favorable on the trailing side. 

Near the tip radius in the ” interference region" near the trailing 
edge at x/c = 0,7 and 0.9, the turbulent intensities are observed to 
increase significantly compared to the free stream values (e.g. 

Figures 65, 69, 73), As explained in the section on Kean Velocity 
Profiles, the two blade boundary layers have completely merged with 
each other and their interaction with the annulus wall produces a large 
turbulent mixing zone in the central region of the passage. Gorton's 
(23) result also confirms this observed behavior for a three bladed 
inducer. Unlike the tt^o-dimensional fully developed boundary layers in 
stationary passages or pipes, the streamwise velocity profile in this 
interference region shows a dip with inward radial flows. There exist 
finite mean velocity gradients even in the raid passage as shown In 
Figures 48 and 49 at R = 0.98 and result in an increased production of 
turbulence in this region. 

5.3.5 Turbulent Shear Stresses 

Figures 53 to 73 also shovj the blade to blade distribution of all 
the three turbulence shear stress correlations derived from the hot wire 
measurements described earlier in Chapter IV. It should be remarked 
that the mean velocity, turbulence intensity and shear stresses are 
measured simultaneously at the same location. Expected error in the 
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measurement of fluctuating voltage (Appendix B) cross-correlations is 
less than 6.7% for ^^^2 ®1®3 9,1% for «2®3 i^®spectlyely. This 

should result in a maximum anticipated error of less than 9% in velocity 
cross-correlations or the turbulent shear stresses. Because of other 
miscellaneous errors the shear stress data should be used with caution. 

The turbulent stresses are small in the center of the channel 
(except in the interference region) and Increase in magnitude on 
approach to the blade surface, and this trend is similar to turbulence 
intensity components. This is brought about by the production of these 
stresses by the mean velocity gradients as shown in transport equations 
given in Appendix C (Equations C-4 and C-6) and is common to all wall 
bounded shear flows. 

However, unlike two-dimensional boundary layers where only the 
streamwlse turbulent shear stress component (-u*v*) is present, the 
data clearly shows that the other two stress con^onents (viz. (-u*w*) 
and (-v’w*)) are also of the same order of magnitude as (-u*v*). In 
three-dimensional stationary boundary layers the component (-v*w*) is 
present due to finite value of the production term in its transport 
equation (Equation C-6) in Appendix C. However, the correlation (-uV) 
which represents a cross-correlation normal to boundary layer is found to 
be negligible because the production term due to mean velocity gradients 
are small in this direction. In the present case, in addition to mean 
velocity gradients normal to the boundary layer, there are large mean 
velocity gradients in radial and chordwlse directions (Figures 45 to 49) 
and result in extra production terms. - — 

It is evident from Equations (C-4), (C-5) and (C-6) of Appendix C, 

that the effect of rotation in our case (where = — QsinB* — Qcos^, 

y * 



3 varying from 78” at the hub to 84” at the tip) Is to Increase the 
streanwlse component (-u^v*) of stress an<j[ decrease the other two com- 
ponents (-v'w' and u'w'). The experimental Investigation on a single 
blade (Ref. 17) as veil as the present set of data for the streamvlse 
component (-u'v*) seem to confirm this trend. It Is also evident from 
Equations (C-4) to (C*-6) that the curvature terms which are small in 
this particular case, seem to have effects opposite to that of rotation 
In the present Investigation. The effects of curvature for a swirling 
jet are discussed by Pratte and Keffer (47). 

As explained In the section on turbulent intensities, only the 
data near the mid radius for x/c > 0.3 In Figures 58 to 65 is free from 
leading edge Incidence and end wall effects. It represents the turbu- 
lence behavior in the turbulent boundary layers on the blades of turbo- 
machinery rotors operating at designed flow conditions. 

It Is seen that the correlations u'v* , v'w' and u'w* follow the 
same trend on the leading side which means that the streamvlse, radial 
and blade normal turbulent shear stresses behave similarly. On the 
trailing side, they follow the opposite trend l.e. (u’v*) and (^w^) 
correlation has values In between the two. This difference Is 
attributed to the effect of chordwlse pressure gradients (Figure 44) in 
the two cases. On the leading side, for x/c > 0.3, the chordwlse pres- 
sure gradients are small and therefore production terms In Eqtiatlons 
(C*-4) and (C-6) due to mean velocity gradients In streamvlse directions 
are small and may be neglected. The major shear production terms are 
V* 3u/9y and v’ 9w/3y In (C-4) and (C-6) respectively and normal 
three-dimensional boundary layer behavior Is observed. On the trailing 
side, the chordwlse pressure gradients are favorable and large. 
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This results in production tertns of opposite sign in two cases resulting 
in the observed trend. 

In the interference region (Figures 58, 65, and 73), the turbulent 
stresses become large due to production by mean velocity gradients. 

The magnitudes of ("uW') and (-u*w*) are once again opposite in sign. 
This is due to the fact that streamwlse velocity gradient (3u/9y) has 
changed sign with Increase in distance (y) normal to the blade surface 
vjhile (9v7/3y) has the same sign because radial flow is Inward. The 
stresses increase almost linearly with y in this region, except where 
the mean velocity gradients go to zero, i.e, at the position, where u 
has a minimum value and inward radial flow shows a maximum. 

’5.3.6 ’ Turbulence Spectra 

Spectra of all the six components of Rejrnolds stresses and 
turbulent kinetic ener^ were measured at two stations at x/c == 0.5, 

R = 0.94 and distances respectively of 3/32 inch (8 = 2,7°) and 5/32 
inch (8 = 4.5°) from the leading side of the blade surface. The main 
objective of the spectral measurements is to find the cut-off frequency 
(to eliminate electrical noise signals If any' above that frequency) for 
the 2-charihel Butterworth filter used' in the hot wire set up in Figure 
11. This cut off frequency is characteristic of the flow as well as the 
limitation of- the hot wire anemometer used in the measurement. It was 
also desired to check the contribution in the hot wire output from any 
spurious electrical signals at line frequency, blade passing frequency 
and natural frequency of vibration of the probe. One of the main 
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objectives of this measureinent is to discern the influence of rotation 
on turbulence structure. 

Measurements (Figures 74 and 75) clearly show that a cut off 
frequency of 5 kHz is quite suitable for turbulence intensity and shear 
stress measurements and the contribution to electrical signal at wave 
numbers higher than 25 cm ^ is negligible in the present flow conditions. 
Also no spurious electrical signal was observed at line frequency, or 
blade passing frequency or natural frequency of the probe below 5 kHz. 

The behavior of turbulent intensity spectra shown in 74 and 75 reveal 
that the inertial subrange where the spectrum follows a (-5/3) law is 
very small and the spectra falls rapidly with increase of wave number at 
large wave numbers. On the other hand, the spectra of turbulent shear 
stresses shown in Figures 74 and 75 reveal a (-5/3) law behavior 
between wave numbers 2 to 20. A detailed spectral measurement both on 
the leading and trailing sides and at various chordwlse locations would 
be required to correlate this behavior with effects of rotation and 
curvature in rotating blade boundary layers encountered in axial flow 
turbomachinery . 
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Figure 74. Spectral Distribution of Turbulent Intensity and Shear Stress Correlations 
at R = 0.94, x/c * 0.5, y"*" 50 (Throttled Case). 
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VI. CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE WORE 
6.1 Conclusions 

The theoretical and experimental Investigation of three 
dimensional rotating turbulent boundary layer In a helical passage of 
a turbomachlnery rotor show substantial differences be twaen their 
stationary counterpart as well as from those on a single rotating 
blade. The rotation Induced Coriolis force markedly changes the 
boundary layer growth, the skin friction, the extent of thret$~ 
dimensionality and the turbulence flow characteristics In a rdtatlng 
channel. It should be pointed out that these detailed quantitative 
measurements on the turbulence and growth characteristics of blade 
boundary layers In a turbomachine rotor have been performed for the 
first time. These Investigations reveal the Inadequacy of the data 
obtained from cascades In calculating or predicting the performance of 
a turbomachine. Following important conclusions can be drawn as a 
result of this Investigation. 

1. A collateral region exists very near to the wall, even for 
highly three*-dlmenslonal rotating turbulent boundary layers with radial 
and streamwlse pressure gradients. • 

2. The law of the wall Is found to be valid up to y^ ^ 200. 

3. A Rower Law representation of the chordwlse velocity components 
fits the data well provided the radial and streamwlse variation of the 
exponent ((H‘>l)/2) In the chordwlse velocity profile Is taken Into 
account. Cole's representation is possible but the wake ftmctlon Is 
different from the sinusoidal form used for the stationary turbulent 
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boundary layers, indicating that the entrainment process in rotating 

f 

boundary layers is affected by Coriolis and centrifugal forces. 

4. Mager’s representation of the outward radial velocity (with 
measured limiting streamline angle) fits the data reasonably well for 
boundary layers near the mid radius but departs considerably when the 
annulus wall boundary layer or the boundary layer from adjoining 
blade interferes. 

5. The velocity profiles in the outer part of the boundary layer 
agree well with the expression, 

= Cl / ^2/k - sinB) sin3/RG - (u/U^)^ 

provided the variation of Ci with x is allowed. In the developing 
region, for open throttle case, where the boundary layer edge velocity 
is f2r, w/U^ = 0.2 J (2-sin3)sinB-(u/U^)^ gives the best fit. It is also 
clear that Johnston straight line fit 

w/U = C_ (1 - u/U ) 
e z e 

is inadequate for rotating axial flow fluid machlneiry, such as axial 
flow compressors and turbines. Equation (31) clearly shows that the 
cross flow is strongly dependent on the stagger angle, local chordwise 
velocity and the radial pressure gradient. 

6. Velocity profiles in the interference region near the tip, 
where the blade and annulus wall boundary layers produce an Intense 
mixing region, cannot be represented by simple similarity forms. Radial 
inward flow and appreciable reduction in the streainwise component of 
velocity characterize this flow. Equation (31) gives the best fit 
except near the center of the channel section. 
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7. There is a large reduction in passage averaged strearawise , 
velocity u downstream of the channel due to increased radial migration 
of the flow caused by the pressure gradients. The decrease is very 
large near the tip in the interference region and is accompanied by 
radial inward flows and reversal of shear s treses. This phenomenon 
plays an important part in producing increased secondary flow losses in 
an actual rotor. Its existence and its quantitative measurements have 
not been reported before. The associated intense mixing and the 
corresponding increase in end wall losses must be accounted for in 
calculating the performance of turbomachinery rotors, 

8. The observed values of momentum thickness lower than 

those for a single rotating blade measured by Lakshminarayana et al, (17) 

due to larger values of the limiting streamline angle (e ) and the 

w 

resulting larger radial velocities. 0^^ is higher on the trailing 
side compared to those on leading side. The predictions based on 
momentum integral solution agree reasonably well. 

9. The observed and the predicted values of limiting streamline 
angle are much higher than the single blade values and agree well except 
near transition region. This is partly due to errors in the location of 
transition and partly due to numerical errors involved in assumed sudden 
transition from laminar to turbulent flow equations. 

10. Variation of shape factor (H) is small in open throttle case 
(1.10-1.3 on the trailing side and 1.25-1,5 on the leading side) except 
near hub and tip where the end wall boundary layers interfere. For a 
fully loaded rotor, H varies appreciably. Its variation can be incorpor- 
ated in the prediction of boundary layer growth by using a modified form 
(Equation 92) of Head’s entrainment equation. 
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11. The skin friction coefficient for the rotating channel is 
found to vairy with three important three-dimensional boundary layer 
parameters ©^i* proposed correlation for the skin 

friction coefficient based on our limited data in rotating channel 
with mild pressure gradients is 


C. = 0.172 R. '^"268 ^q- 0.678H ^4.3^52 /e (x-x^/C ) 

f 0^^ w t 

The term in the bracket is due to three-dimensional effects. Variation 


with Rq and H is found to be the same as that for a stationary channel, 

o 


11 


12. Turbulence intensity measurements show that the radial 

— T 

component (w’ ) is larger than the strearawise component (u* ) in a 
rotating turbulent boundary layer. This trend is opposite to that in a 
stationary channel and is due to rotation. 

13. Turbulent shear stress measurements show that in three- 
dimensional rotating turbulent boundary layers, all the three velocity 
cross-correlations are of the same order of magnitude inside the boundary 
layer. They behave similarly if the pressure gradients are negiglble 
and can have opposite signs depending on the nature of chordwise and 
radial pressure gradients. 

14. Both the turbulent intensity and shear stress measurements 
in the interference region (near the tip) show a large increase in 
magnitude. The turbulent stresses vary almost linearly in the direction 
normal to boundary layer and have opposite signs compared to the values 
near the blade surface. This large gradient in turbulent intensity and 
shear stresss is due to an Intense mixing zone due to interaction of the 
two boundary layers resulting in increased flow losses in this region. 


172 



6. 2 Recommendations for Future Work 


The present study on the three-dimensional turbulent boundai^ 
layers in the rotating blade passages of the four b laded axial flow 
Inducer has shown that there is an urgent need to do both analytical 
and experimental work in other rotating fluid machinery, such as 
compressors, turbines and propellers. The turbulent characteristics, 
the skin friction and the boundary layer growth are significantly 
affected by both rotation and curvature and quantitative measurements 
of these effects in rotating frame of reference is essential to predict 
the losses and the aerodynamic performance of a turbomachinery. As a 
result of this study, following recommendations are proposed for further 
research work in the area of rotating turbulent blade boundary layers: 

1. A complete solution of mean velocity profiles and the major 
turbulent stresses in the blade boundary layer using Bradshaw’s (25) 
three-dimensional boundary layer program for stationary flows, modified 
to include the effects of rotation and curvature as outlined by Anand 
et al. (20). 

2. A prediction method for the boundary layer growth parameters 

(0--, £ , and H) by matching analytical solutions for cross-flow 
11 w 

velocity profiles in the inner and outer regions of a turbulent 
boundary layer, and using the momentum integral solution technique 
described in Chapter II for arbitrary pressure gradients in a turbo- 
machinery rotor. The entrainment Equation (92) should be modified to 
Include the effects of free stream turbulence as discussed by Evans and 
Horlock (48) and the effect of cross-flow on the Head's (32) entrainment 
coefficient. 



3. A conq) rehens Ive study of Interference effects of annulus wall 
on the blade boundary layer characteristics both In the stator and the 
rotor blade row of a turbomachinery. 

4. Accurate boundary layer measurements close to the blade surface 
to verify the extent of collateral region In the Inner layer as 
suggested by Equations (24) and (25) . 

5. Use of boundary layer probes with high spatial resolution, 
especially In the measurement of turbulent stresses and the cross- 
spectra. 
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APP^DIX A 

Specification p£. Hotwire Probe and AneTOoeter and 
Turbulent Scales of Flow 

A.l Specif iCatibtis of the Hotwire Probe ■ 

All the hotwire measurements reported In Chapter V have been 
carried out using a miniature triple sensor hotwire probe manufactured 
by Thermosys'tems • Inc. The configuration of the probe and its sensors 
is shown in Figure 14. The probe Is specially suited to boundary layer 
measurements. The probe body Is 0.08 inches thick stainless steel tube • 
eleven inches long and elliptic In section with major and minor axes 
0.25 inches and 0.125 inches respectively. 

The three hotwire sensors of^ the probe are supported by six gold 
plated stainless steel tapered needles, 0.008 Inch tip diameter about 
0.35 inch away from the probe body facing the direction of the flow, so 
that the Interferenca-ef fectSi due to probe body on sensors are minimized 
The sensors are made of copper plated tungsten wire, 0.00015 Inch In 
diameter* The total- length of each Sensor Is 0.0624 Inch but the etched 
length varies from 0,025 to 0,040 Inch. The properties of hotwire 
sensors are tabulated in Table 1, their direction cosines In Table II 
and the influence coefficients In Table III, below: 




Table 1 


Important Properties of Hotwire :^Seti8or Meterlal 


b . 

Property 

Symbol 

Value 

Unit 

1 

Temperature Coefficient of 
Electrical Resistivity 

,(b) 

5.2 X 3.0^^ 

/O 

>: ;v. -‘./c 

2 

Density 


19.0 

gm/cm^ 

3 

Specific Heat 

c 

w , , . , 

0.032 

BTU/lb. 

A , 

Thermal Conductivity 

,k 
■ w 

0.A76, 

gm calorie/ 
sec °c 


• Table II 

Direction Cosines of Sensors with Reference Coordinate System of Flj^ure 



Hotwire 1 

Hotwire 2 

Hotwire 3 


Angle 

Angle In 

Direction 

Angle in 

Direction 

angle in 

Direction 


Orientation 

degrees 

cosine 

degrees 

cosine 

degrees 

cosine 


Ax - X 

135° 

-0.707 

118.5 

-0i477 

232.5 

-0.609 ’ 


Ax - y 

A5° 

0.707 

30.5 

0.862, 

141.9 • 

-0.787 


Ax “ r 

45° 

0.707 

80.0 

0^74 

85.0 

.087 



Ax - X, Ax ~ y, and Ax ' r represent Included angles between the wire 
axis (Ax) and x, y, and r axes» In the planes foi^d by wire, axis and 
X, y, r axes respectively* 
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Table III 


Influence Coefficients of Sensors Used In Hotwire Equation 

(99) In Chapter IV for Gorton’s (23) Computer Program 


Hotwire 

“u 

^£2 

“£3 

1 

0.720 

0.0 

0.670 

2 

0.882 

0.455 

0.092 

3 

0.807 

-0.558 

0.062 


A. 2 Specifications of 2-Dual Channel Constant 

Temperature Hotwire Anemometers 

The anemometry circuitry has been designed to process large 
nonlinear fluctuating signals by Wyngaard and Lumley (42) with a flat 
frequency response up to 20 kHz. The electronic circuit diagram is 
shown in Figure 76 and the details of conq)onents are given in Table IV. 
Since the spectral measurements show that the magnitude of output 
signal reaches the noise level (0.02 mv - anemometer in operation for at 
least one hour) for Oi < 5 kHz in the actual experiment, this anemometer 
is found suitable for all the hotwire measurements reported in Chapter V. 
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Table IV 


Anemometer Components . ~ 


Resistor No. 

Value 

Q 

Watt 

Composition 

Number 
required 
for two- 
channel 


lobfi pot 

1 turn 


2 

R3, R3 

100 

1-min 

wire wound 

4 


5000 

1/2 

metal film 

2 


2000 

1-min 

wire would 

2 


1000 

1-mln 

metal film 

2 


lOK 

1/2 

metal film 

3 

‘‘10 

50 K 

1/2 

metal film 

3 


1000 

1/2 

metal film 

3 

«12 

2000 

1/2 

metal film 

3 

■^3 

500 

1/2 

metal film 

2 

^14 10 

K Potentiometer 

10 turn 


2 


lOOK 

1/4 

carbon 

2 


Capacitor 

Value 

Quantity 


1 yf 

4 

S 

560 yf 

2 

S 

.002 yf 

2 


2.2 yf 

3 

s 

68 yf 

3 

s 

1500 yf 

2 

filter capacitor 6.8 yf- 35 V 

16 

filter capacitor .01 yf-low volt. 16 
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A. 3 Scales of the Flcjw 


The relevant range of flow parameters for estimating various scales 
of turbulence for the present flow study are: 

1. Free stream velocity In a relative frame of reference ■ 30 to 
60 ft/sec, 

2. Boundary layer thickness = 0,04 ft at x/c » 0,1 to 0.12 ft at 
x/c > 0,5. 

3. Skin friction velocity = 2 to 3 ft /sec, 

4. Overall turbulence level - 5% to 15% Inside the blade boundary 
layers . 

5. Reynolds number based on tip radius and peripheral velocity at 
the tip = 6.67 X 10^, 

An estimation of turbulent length scales Is very Important from the 
point of view of categorizing the type of turbulent flow, understanding 
Its basic structure and In calculating the various errors Involved In the 
measurement of turbulent flow correlations, which are described In 
Appendix B. The three In^ortant scales of turbulence are: 

1. Large eddy length scale (£) 

2. Taylor’s microscale (^g) 

3. Kolmogorov scale (ri) 

Following Tennekes and Lumley (49) , the large eddy length scale Is 
proportional to the mean flow strain rate. It Is of the order of the 
mixing length In the case of turbulent boundary layer flows. H Is 
therefore ''^0.2 Inch at x/c = 0.1 and 0.6 Inch at x/c > 0.5. 

The Reynolds number (u’Jt/v) Is an Important parameter In determining 
the structure of turbulent flow. For the present case Its value lies 
between 500 at x/c a 0.5 to 1500 at x/c >0.5. Here u* represents 
the turbulent Intensity. 
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Taylor's microscale (1 ) is defined as: 

g 

e = 15 Vu'^/A ^ 


(A-1) 


where e is the turbulent dissipation rate (£). An estimate of e 
according to Tennekes and Lumley (50) is obtained by 


e = n'^/l (A-2) 

2 3 2 3 

Therefore e is 1000 ft /sec to 1500 ft /sec and X is 'v 0.04 

g 

inch to 0.07 inch at x/c = 0.1 and x/c > 0,5 respectively. 

Kolmogorov scale (0) which represents the small scale eddies in the 
dissipation range is defined as 

n = (A- 3) 

-3 

The Kolmogorov scale is 3 - 4 x 10 inch. 

In the spectral representation of Reynolds stress tensor components, 

the concept of wave number is very useful. The wave number is defined as 

K = 2'n-w/u - 1/il (A-4) 

e 

where is the frequency at which the contribution due to a correlation 
between two fluctuating velocity components to the Reynolds stress is 
being considered, u is the mean flow velocity at which an eddy of 
length corresponding to that wave number is convected. 

The wave number is usually represented in centimeter units. The 
wave number corresponding to an eddy of wire length is 13 cm ^ and the 
wave numbers corresponding to an eddy equal to distance between wires 1 
and 2 or wires 1 and 3 (Figure 14) is 10 cm The wave number 
corresponding to distance between wires 2 and 3 is 6.4 cm These 
wave numbers are required in calculating errors in Reynolds stress 
correlations described in Appendix B. 
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APPENDIX B 


Experimental Error Analysis and Accuracy of Data 

The various sources of error in the measurements and methods of 
incorporating them in the data are described in this Appendix. 

B . 1 Limiting Streamline Angles 

The following are the sources of errors in the measurement of 
limiting streamline angles by the ammonia trace technique, 

1. Finite velocity of the injected ammonia and extent of collateral 
region. 

2. Finite dimensions of the static pressure tap through which the 
ammonia is injected. 

3. Misalignment of the ozalid paper strip with reference direction. 

4. Inaccuracy in measuring angles from the ammonia trace due to 
increased blurredness in the trace with distance due to turbulent 
diffusion effects. 

Items 1 to 3 are controllable and are minimized by experience and 
proper design of the ammonia transfer device. However item 4 is the 
main source of inaccuracy. An error in angles up to ± 2° is possible. 
Therefore the percentage error in the small cross flow case (where the 
limiting streamline angle is ^ 10°) can be as high as 20%. 

B, 2 Skin Friction Coefficient 

The Preston tube measurements for the blade surface shear stress are 
liable to the following errors: 
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1. The Inner diameter of the probe should be small compared to the 
exterior law of the wall region. Near the leading edge the * law of 
the wall* region Is very small and results In some error. 

Misalignment of probe compared to the actual flow direction at 

the wall . 

3. Measurement of effective radius of the probe location from the 
rotor axis, when the probe Is Inclined. 

Error due to Item 1 Is found negligible, as demonstrated by the 
measurements with different tube diameters. The error due to item 2 is 
minimized by using a rounded tip Preston tube. The last error due to 
Inclination of Preston tube depends upon the magnitude of limiting 
streamline angle and the radius of measurement. It Is maximum at the 
hub for large pressure gradient cases. For a one Inch tube, tilted 
up to 45^, the total anticipated error In measurement at hub Is 6%. 

B. 3 Mean Velocity Profiles 

The sources of error and experimental accuracy depend both on the 
type and geometry of the probe and Its sensors. Both the three hole 
disk probe and the three sensor hot wire probe, have been used for mean 
velocity measurements. Probable sources of error from disk-type probe 
measurements are: 

1. Nonuniformity in the flow. 

2. Proximity of the probe to the blade surface. 

3. Finite dimensions of the sensor area. 

4. Inclination of the streamline to the total pressure hole. 

5. Misalignment In reference direction (x) , which Is parallel to 
the blade surface. 
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All errors except (5) are due to Inherent characteristics of the 
flow or the probe. They are minimized by Incorporating the following 
corrections to the data: 

1. Correction due to nonuniform plane shear flow In the radial 
direction (Figure 77): Since the calibration of the probe is 

conducted in a wind tunnel where the flow is uniform, an error due. to 


shear flow arises, when the probe scans through the blade boxmdary layer. 
The magnitude of this error in angle measurements depends upon the 
stagnation pressure gradient in the relative flow and the distance 
between the two pressure holes used for the flow angle measurements. 
Error in angle measurements (e^) due to shear flow in radial direction 
at any measuring point can be proved to be 


k i = k- ^ 

1 8r u 1 r 


(B-1) 


where Z is the radial distance between the two holes and is a constant, 
which is found to be nearly equal to 1.0. For the present case, Z « 0.175 
inch, u ^r and r varies from 9 inches at hub to 18.25 inches at tip. 
Hence the error in angle measurement (e^) 1.0*^ at hub and 0.5° at tip. 

The error due to finite dimensions of the sensor area is directly 
related to the distance between the holes. From the practical point of 
view, for 1/32 inch diameter holes, a minimum distance of 5/32 inch is 


required and the design value of 0.175 inch is quite appropriate. 
Maximum anticipated error in the velocity profile measurements 'v* 1% at 
the h:ib. It is partly corrected by using Equation (B-1). 

The second source of error due to shear flow is the displacement 
of total stagnation pressure point from its true position. LighthllL's 
(50) expression for the displacement (A) of the effective center of 
pressure from the probe axis for a cylindrical probe of diameter d is 
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(B-2) 


is. d 

I “ 0*^5 C - 3.05 5^. where ? - |^ (|;j) 

Here the shear gradient In the direction (y) normal to the boundary 
■ . 3u 3u 

layer Is used, because . In the present case, d ■ 0.031 Inch, 

C ~ 0.02, 0.009. The traversing step In boundary layer measurement 

Is 0.06 Inch. The error due to this shift Is therefore at the most 1.8X. 

2. Correction due to proximity of the blade surface: The effect 

of the blade wall on probe measurements Is to displace the effective 
center of measurement from Its true position In the direction of 
Increasing distance from the blade. According to Thwaltes (51) the 
maximum displacement (A) due to this effect Is given by ~ ^ where t 

Is the thickness of the probe. The disk probe Is 0.03 Inch thick. 
Therefore the maximum displacement Is A < 0.015 Inch. Even though It 
Is small. Its affect must be accounted for the measurement closest to 
the wall. This correction has been Incorporated In the final analysis 
of the data. 

4. Error due to Inclination of the streamline to the total pressure 
hole: Total pressure measurements will be In error If the streamline Is 

Inclined to the total pressure hole of the disk probe. This error Is 
eliminated by using a calibration curve for the total pressure measure- 
ments. This Is carried out by traversing the probe Inside the potential 
core of a 1 foot diameter open jet. The total pressure Is measured for 
different angular positions with respect to jet axis. The normallaed 
difference between true and observed total pressure Is plotted against 
flow angle for different values of flow velocities In the tunnel. The 
single correction curve obtained by this procedure (Figure 78) is used 
to correct the total pressure measurements as follows: 






1 + 




<B-3) 


om 

where h and h are, respectively, the measured and actual values of 
om o c 

total pressure head of air In Inches of water and Is the correction 

from the calibration curve at any given angle a (Figure 78) , It Is to 

be noted that the flow angle Is also to be corrected because the 

normalizing dynamic head In the angle calibration curve (Figure 13) uses 

the uncorrected measured head (h ) to calculate the dynamic head. Thus 

om 

both OL and h values are corrected by Iterative procedure. It may be 
om 

remarked that the error In total pressure for a flow angle < 10^, Is 
less than 3% and hence the correction may be neglected. 


B.4 Errors In Hotwire Measurements 

In the case of measurements with the three-sensor hotwire probe, 
the following sources of errors may occur. 

1. Inclination of the wire to the flow streamline (deviation from 
cosine law) . 

2. Geometry of prongs and probe body (Inviscid flow and heat 
transfer effects) , 

3. Finite distances between wires (spatial resolution of the probe) . 

4. Finite dimensions of Individual wires (^/d, end effect, spatial 
resolution, etc.). 

5. Temporal resolution or thermal Inertia of the wires. 

6. Aging, oxidation and contamination of the wires. 

7. Ambient temperature drift. 

8. Proximity of wire to the wall. 
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9. Probe body vibration due to rotation and flow induced excitation. 

10. Spurious signals from the power line (60 Hz) and blade passing 
frequency. 

11. Finite sampling time. 

12. Measurement of wire angles with respect to the reference 
coordinate system. 

13. Misalignment of the probe with respect to the reference axis. 
Errors 1, 7, and 8 are eliminated by using correction factors from 

available hotwire literature (52, 53, 54) and calibration in a wind 
tunnel. The probe and sensor errors 2 to 6 are minimized by proper 
selection of the probe dimensions, support needles and the sensors. 

Error 9 is eliminated by proper clamping, error 10 is minimized by 
reducing the number of ground loops to a minimum. Error 11 is Inherent 
and can only be reduced by increasing the sampling time. Errors 12 and 
13 depend on the accuracy of angle measuring equipment and aligning 
devices respectively. 

Estimation and minimization of all the above mentioned errors and 

methods of correcting the data are described below. 

1. Inclination of wire to the flow streamline: Inclination of a 

hotwire to the flow streamline causes a change in effective cooling 

velocity on the wire and is expressed by Hinze (52) as, 

2 2 2 2 

U .. = a/ + U (B-4) 

eff N a 

where and are components of total flow velocity vector respectively 
in normal and axial directions of the wire respectively. K is a 
correction factor, which according to Friehe and Schwarz (55) depends 
mainly on the il/d ratio of the hotwire. 
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2. Geometry of prongs and probe body; The prongs and the probe 
body affect the hotwire measurements in two ways. Firstly, they change 
the flow velocity vector near the wire due to obstruction of the flow 
(inviscid effect) and secondly it affects heat transfer rate. The major 
design parameters which influence this error are prong to wire diameter 
ratio, prong length to diameter ratio and spacing between the prongs. 
From Appendix A, the values of these are respectively 20, 40 to 50 and 
2ram. Following DISA (54) , the maximum anticipated error for 45° 
inclination of streamline from the calibrated position (a =» 0) is 1.5%. 
Since the flow vector inclination from calibrated position in the 
measurement is less than 45°, this error is negligible. 

3. Spatial resolution of the probe: The finite distance between 

the wires or the total sensing area, can be a large source of error, 
particularly in nonuniform three-dimensional flows. Since the hotwire 
measurements by individual wires are made at 40 points from one blade 
to the other, it is easy to interpolate the measurement at a single 
point from the observed measurements at adjacent data points. Since 
wire 1 in Figure 14 coincides with one of the reference axis (x) the 
measurements of the other two wires are interpolated to this point. 

Since the distances between wires is small (maximum value is 0.060 inch) 
as compared to measured boundary layer thickness (y 0.5 inch), a linear 
interpolation scheme is used. Unfortunately this correction cannot be 
applied to the measurements of cross-correlations and the corresponding 
error in their measurements is derived as follows: 

The error in a cross-correlation due to spatial resolution is 
nothing but the minimum eddy size the probe can correlate because of 
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separation of the two wires, 
these distances are ° *^^.3 
and described In Appendix B. 
Lumley (56) as 


The wave numbers (tc) corresponding to 

-1 -1 
= 10.0 cm and * 6.4 cm as defined 

The error e In the correction Is given by 


I E^^(k) dK 


‘Ij 


e(u,»uj«) 


J E^j(k) dK 


- u^*u^*(k) 


(B-5) 


o 

The turbulence can be taken to be locally isotropic for K 5 cm An 

m w 

estimate of correlation and hence the error may be obtained by assuming 
that the spectrum of E follows the following expressions as discussed 
by Wyngaard (57), 


Eij(K) * constant C^(0 < k < 0.1) 

* constant C^/k (0.1 < K < 5,0) 

2/3 -5/3 

= constant e k (5.0<k<«>) 


(B-6) 


The correlations u^'u^* are directly related to measured fluctuating 


voltage correlation described in Chapter IV. Using these 

relations, the maximum anticipated errors in the measured voltage 
correlations are: 


Error in ®1^3 ®°^’^®l®tions = 6.7% 


Error in ©2^3 co^^relation 


= 9.1% 


So the anticipated maximum error in cross-correlations due to spatial 
distance is ^ 9%. 

4. Spatial resolution of the wire: The finite dimensions of the 

hotwire cause the following errors: 

a. Error in heat transfer rate parallel to axis of the wire 
(il/d) ratio) . 
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b. Error in heat transfer rate due to cold length effect 

c. Error in measurements due to spatial resolution In 
turbulence data. 

d. Error In measurement of mean velocities In nonuniform flow. 

e. Thermal Inertia error In fluctuation levels at different 
frequencies . 

The first error Investigated by Frlehe and Schwarz (55) can be 
Incorporated as outlined in Equation (B-4) by making K as a function of 
jl/d ratio, where ^ and d are respectively the etched length and diameter 
of the wire. The H/d ratio varies from 133-300, Following Frlehe and 
Schwarz (55), K ranges from 0.25 to 0.15 respectively. 

The cold length effect or end effect Is due to transfer of heat to 

prongs at the end of the wire. Following Hlnze (53), the cold length' 

is given by 

where e Is a conversion constant between heat and electrical units. For 
c 

heat in calories and electrical energy in joules, its value is 4.2 fc^ 

Is the heat conductivity of the wire material, b Is the temperature 
coefficient of the electrical resistivity of the wire. is the 
electrical resistance of the wire at 0*^C. A and B are hotwire constants. 
I is the value of electric current flowing In the wire to maintain a 
constant temperature, when the fluid flow velocity is U, For A * 0.35, 

A + B » E = 0.45, I = 80 m A and using the specifications of the 

wire from Appendix B, the cold length 5.^ is given by 

Z I 

-| = 60 and -^ = 0.25 to 0.45 
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This error Is eliminated by calibration' of each wire over the complete 

flow range and using the local values of slope B In the King's law. 

The error due to spatial resolution of wire is discussed by 

Wyngaard (57). Since the Kolmogorov length scale ri varies from 0.003 

inch to 0.005 Inch. The corresponding wire length to Kolmogorov length 

scale ratio (Vo) Is 6. This gives about 10% error in intensity 

spectrum at a wave number given by ^^^5, = 1,2 or “ 10 cm ^ where 

K. Is the wave number where the measured spectrum of mean turbulent 
im 

energy falls by 10%. Since the magnitude of spectrum falls to very 

small values for K>10 cm and area under spectrum for k . ^10 Is about 

Im 

90%, the total error in turbulence Intensities Is < 1%. 

/ 

Since the sensor wires are in a nonuniform three-dimensional flow 
field the effective cooling velocity varies over the entire length due 
to their finite length. The calibration curve, on the other hand Is 
based on uniform velocity field over the wire. The measured output 
voltage therefore doesn't correspond to mean velocity at the center of 
the wire. Similarly an error is also made In turbulence measurements. 

The error in turbulence measurements due to nonuniform flow Is given by 
Hlnze (52) as 

e = C (e ) measured (B-8) 

~~2 ' 
where e is the corrected voltage and C is the correction factor given by 

-1 

C ^ (1 - ^^) ^ (for l/\ < 1) (B-9) 

6X 2 8 

g 

where Z Is the sensor length and X is the Taylor's micro-scale of 

g 

turbulence. For the present case, I ^ 0.030 Inch and X » 0.040 to 0.068 

g 

inch. Therefore the correction factor to rms voltages varies from 1.03 
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to 1.09. This factor If neglected Introduces a maxlmuin error of 4.5% In 
the Intensity measurements. 

5. Thermal Inertial of the wire: The temporal error due to 

thermal Inertia at any frequency U) Is given by Hlnze (52) as, 

e(oj) *• 1/ 1 ^ + 0) where M Is the time constant of the 


wire 


M 


e C e C (9 -0J 

c w w f c w w f 


b I^R R- 
o f 


I^R. 


(B-10) 


where C Is the heat capacity of the entire wire, R and R^ are the 
w w r 

values of wire electrical resistance at temperatures 6^ and 6^ 
respectively. 6^ Is the temperature to which the wire is heated and 6^ 
Is the fluid temperature. 

For tungsten wires used In the experiments with over heat ratio of 

-4 

0.8, and using the specification data from Appendix A, 1.4 x 10 
secs. Therefore Thermal Inertia error In measured fluctuation levels 
for 1 kHz and 5 kHz are. 


e (1 kHz) » 0.98% 
e <5 kHz) - 25% 

6. Aging, oxidation and contamination of the wire: The errors due 

to aging, oxidation and contamination depend on the total operation time, 
over heat ratio and the flow environments respectively. These errors are 
minimized by calibrating the probe before and after the experiment, and 
Incorporating the changes In calibration constants by Interpolating 
according to the fraction of total time used In each measurement. 

7. Ambient temperature drift: Change In ambient temperature causes 

changes Ih the heat transfer rate of sensors resulting In corresponding 
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change In measured output voltages from the anemometer. This effect can 
be corrected In the following manner. 

2 r** 

The constants In the hot wire equation E = A + B are given 

m 

by Hinze (52) as, 

A = A k^(Pr -0.) (B-11) 

o r m w I 

and 

B = B k Re,°*^(Pr )^*^^(0 - 0J 

o f d m w f 

where 0^ is the wire temperature when heated, subscript m refers to 

values at mean temperature equal to (6 + 0.)/2, and 0- is fluid 

w r X 

temperature, is thermal conductivity of wire at fluid temperature. 

and B^ are functions of wire dimensions and those properties of the 

wire that are independent of fluid temperature. For common range of 

variation of fluid temperatures of ±10°F, Prandtl number (Pr) is also a 

constant. The thermal conductivity (k ) increases linearly with 

m 

temperature as k = a + a, 0 but k is almost independent of 

m 1 m m a 

temperature. Following Hinze (52) , the effect of ambient temperature on 
two coefficients is given by 


A = A.(l + a6J(0 - ej 

1 r w r 


(B-12) 


B = (9^ - 0p 


where A^ and are constants, a = a^/2 (a + a^/2 0^) *= constant for a 
constant temperature hotwire anemometer. The effect of ambient 
temperature 0^ can therefore be written 


as 


a 


A0, 


A A^ (1 + 

o f o 

A0 

» = Sf 

o w f o 


•)(1 - 


0 -9, 
w fo 


(B-13) 
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where A0_ = 9^ - 0^ Is a small quantity ±10°F, subscript o refers to 
r t ro 

zero flow conditions. For an over heat ratio of 0,8, 9^ = 533°R, a = 
0.0024, 0t/(l+0t8^^) = 0.001 and (1/(6^ - ^ 0,0024 respectively. 

The coefficients A and B are 

A = A^^Cl + 0.001 Aepd - 0.0024 A0^ = A^^(l - 0.0014 ASp (B-14) 

B = B- (1 - 0.0024 A0_) (B-15) 

to t 

The correction factors given by Equations (B-14) and (B-15) have been 
incorporated in the computer program for processing hotwire data. 

8. Proximity of the wire to a wall: The proximity of a hotwire to 

the blade surface causes an error in measurements by increasing the heat 
transfer rate. The increase in heat transfer gives effective cooling 
velocity higher than those in the absence of wall. According to Oka 
and Kostic (58) , the dimensionless difference between apparent and true 
velocity has the form 

^ = f(u^y/v) . (B-16) 

The influence of a wall has been shown by them to be less than (O.lu^) 
for y = u^y/v >6. Since the nearest measurement is taken at y >50, 
error due to this effect is negligibly small. 

9. Probe body vibrations: The vibrations of probe body due to 

rotation or flow induced excitation, if present, can add significant 
amount of spurious signal to the measurements. Even though the probe 
was securely clamped and no visible or other indications of vibrations 
were found (e.g., the wires would break if the probe is vibrating), the 
spectra of the correlation tensor are analyzed at natural frequency and 
its multiple harmonics (lower than 5 kHz) at few representative cases 
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(near the tip radius where the natural frequency of the probe is 
lowest 0(1 kHz)). The spectral measurements reported in Chapter V 
clearly show that there is no spurious signal due to vibrations of the 
probe. 

The natural frequencies of the probe body at various radii are 
obtained by assuming the probe to be supported at one end like a 
uniformly loaded cantilever of elliptic cross-section (Figure 79) , with 
a weight equivalent to the portion of the probe bent at right angle to 
the probe length and supporting the sensors. The effect of probe 
rotation is also included. The natural frequencies in bending mode are 
given by Hartog (59) as, 

u ^ = (U ^)o_„+ a (B-17) 

Ti n 

where (w ^ = 3.03 EI/5, ^ (Hf0,023m) 

E = Elastic Modulus = 30 x 10 psi 

m = mass per unit length of the probe = 0.0018 

M = Point load = 0,00225 lb mass 

—A 4 

I = Moment of Area = 1.7A x 10 in 
a is a constant. Its value according to Hartog (59) is 1.5 for first 
mode of vibration. The natural frequencies based on Equation (B-17) at 
various radii are given by, 

radius (inches) = 18.0, 17.0, 15.75, 14.0 
0 )^ (kHz) = 0.835, 1.4, 1.75, 5.8 

Therefore, for radii lower than 14 inches, natural frequency is greater 
than the cut off frequency (5 kHz) and so does not affect the measure- 
ments. Spectral measurements at tip reveal no spurious signal due to 
vibrations of the probe, if any. 
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Elliptic Section 


M (Equivalent Mass) 



Cantilever Beam of Uniform Loading (m) 


Figure 79. Cantilever Model for Calculation of Natural Frequency 
of Vibration of Hotwire Probe. 
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10. Stray signals from power line and blade passing. frequency: 

Stray signals at line frequency (60 Hz) could be present in the output, 
if there are too many ground loops in the electronic circuitry. The 
number of ground loops is reduced to a minimum and the spectrum measure- 
ments are carried out at 60 Hz. It is found that no significant error 
is present diie to this source. The spectrum is also measured at blade 
passing frequency (190 Hz) and no spike is noticed in the spectrum. 

11. Finite sampling time: The error in hotwire measurements due to 

finite sampling time is discussed by Tennekes and Lumley (49) . The root 
mean square error in true mean value U is given by, 


rms 


(U.J - U) 


- J 2u*^ T, 


/T 


(B-18) 


where T is the averaging time and T is the integral time scale defined 

L 


by 




P(t) dt 


(B-19) 


where p(t) is the correlation coefficient between two values of U 
separated by a time t. According to Tennekes and Lumley (49) the 
integral time scale is of the order of large eddy turnover time, 


T_ 0.4 6/U = 2.7 to 8.0 x 10 ^ secs. 

E 

Also from Hinze (52), for stationary processes, an estimate of the 
integral time scale is given by 


T„ = 2A /U ~ 0.7 to 1.75 X 10 ^ seconds 
E g 

where A is Taylor’s integral scale. The maximum error in the measurement 

s 

of mean velocities for an averaging time (T) is, 

% Error =(100u’ /U)/ 2 t /T 'v 10 x /2 x 8 x 10“^/T 
- 0.40 //t~= 0.07%, 0.06%, 0.05% respectively for T = 30, 
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45» 60 secs. Similarly the root mean square error for turbulent 
components which are almost guasslan with zero mean Is given by Margolin 

(43) as, ^ 

Z(e)^ - 100/ -^ < lOo/^ ^^ T - <J^ (B-20) 

From Equation (B-20), the error in turbulence intensity is IX, 0.85X, 
and 0.7X respectively for T « 30, 34, and 60 seconds. 

12. Measurement of wire angles with respect to reference coordinate 

system: An optical projector of Applied Research Laboratory, The 

Pennsylvania State University, which has a least count In angle measure- 
ments of 0.08° (5*) was used. For angle measurements in the range 45° - 
60° the anticipated error in velocity measurement is 

e =* Aa tana = 0.14% - 0.25% 

13. Misalignment of the probe with respect to the reference axis (x) : 

The probe misalignment with the reference axis. Introduces a significant 

error in the blade normal velocity component v. The measured value of v 

m 

is related to true value by 

V = V ± u sin 6a (B-21) 

m true 

where 6a is the misalignment angle. The maximum anticipated misalignment 
is ±4°. This can therefore add or subtract an error equal to. 6.6% of the 
streamwlse velocity component u. Since the velocity component v is very 
small and goes to zero at the blade surface, the error due to probe 
misalignment is minlmuzed by assuming that v = 0 at the first data point 
nearest to the blade and the corresponding value of misalignment angle is 
used to correct the measurements. 
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APPENDIX C 

Equations of Reynolds Stresses in Rotating Coordinate System 


The transport equations of Reynolds stresses In rotating coordinate 

frame of reference (x, y, r) are obtained from Nash and Patel (31) by 

2 

retaining the principal curvature term ■ sin 3/r and Including the 
additional terms due to rotation. The equations are: 


,.2 


.2 


\ m + +w— +k w) (■= — ^ 

Dt W •' '3t “ 3x 8y 3r xr”’' ^ 2 ’ 


- f (-u>v')+(-u>W)(2By + + 1^ 


II 


xr 

IV V III 


u’ 


— ( 

8y ^ 2 


+ ^ (1^) + V u*V^u' 


VI 


p '3x 
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VIII 


(C-1) 
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p /W’\ _ ,9 I 3 . 3 . 3 x,w’\ 


(-2k^^u - 2fty+ 1^ (.„•„• )+(+2S^+ ^) (17*7^). |£ „.2 


9 /V*w' N , P* 3w* . 2 

- (— 2— >+ ^ ^ '’ ^ w 


(C-3) 


where » -ft sin B, ^ cos B and ft^ - 0 for the flat plate Inducer. 

It should be noticed that the turbulent energy equation (sum of Equations 
C-1* C-2, and C-3) would contain no terms due to rotation. This is due 
to the fact that the coriolis force has no component in the chordwise 
direction. The turbulent stress equations are» 

dT (|?+"lj+’'l 7 +w|^+ W) c-u’v') 

I II I V IV III IV 


(u’V^v* + vW^u‘)+ (u*v*^ + £l (|~+ —) 

P ' p Bx ^ 


(C-4) 


I 



(-U-W) . (|^ + u + V + w I; + k^^w)C-u’w-) 


(-2*'xr“" |j) “•' + w’' (l^xr" + 


V II 


V III 


_[|£ (_uV)+ |H (-v>w’)+(|^+ ||) (-u’w’)]+|^ (u’v'w’) 


II III 


- (-r + ) - V(u V w" + W V U ) 

p 9r 


(C-5) 


5F (-Vw')=(|j-+ u V |^+ w I^K-v'w’) 

- v'2 |a +(2k^^u- 1^ +2fly) (-u’v')-(|^ + |^)(-v’w’) 


V II IV 


III I 


, ,2 ,2s . 9 / i2 , . W*p\ 

•2U (w* -v' )+ *:r-(v" w'+ ■ ) 

X 9y p 


j2Ll (iHl + izl) 

p ^9y ^ 9r ^ 


- V (v'V^w* + w'V^v') 


(C-6) 


These equations are valid for both compressor and turbine rotors. For a 
compressor, Is +ve, 9,^ Is -ve and for a turbine rotor, Is -ve and 


The term on the left hand side Is the total convective derivative 
of a Reynolds stress, terms I to III on right hand side represent 
production of Reynolds stress by mean velocity gradients, term IV Is due 
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to rotation, term V Is due to curvature, term VI represents diffusion of 
Reynolds stress due to gradients of pressure velocity and triple velocity 
correlations, term VII is the pressure strain or equalizing term and the 
last term VIII represents vlsco^ diffusion and dissipation of the 
Reynolds stress. The qualitative effects of rotation terms in the 
transport of a Reynolds stress depends on the directions of rotation and 
are tabulated below for decelerating relative flow in compressors and 
pumps and accelerating relative flow in turbines. 

Table V 

Effect of Rotation 


Transport Equation Rotation Term Effect 


Term Compressor and Pump Turbine 


C-1 

CM 

20 . C-u'w') 

y 

Suppression 

Production 

C-2 

v'^/2 

w’^/2 

0 

No Effect 

No Effect 

C-3 

-20 (“u’w*) 

y 

Production 

Suppression 

C-4 

-u* v’ 

-20 (-w*v’) 

y 

Production 

Suppression 

(if (w*^ - u'^) 
Production 

Suppression 

Production 

is positive) 

Suppression 

C-5 

-u’w* 

+20 (w*^- u’^) 

y 




(if (w’^ - u'^ 

is negative) 

C-6 

-v'w' 

+2fi (-u'v*) 

y 

Suppression 

Production 
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